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Editorial 


POLISH MARITIME RESEARCH is a scientific journal of worldwide circulation. The journal appears as 
a quarterly four times a year. The first issue of it was published in September 1994. Its main aim is to 
present original, innovative scientific ideas and Research & Development achievements in the field of : 


Engineering, Computing & Technology, Mechanical Engineering, 


which could find applications in the broad domain of maritime economy. Hence there are published papers 
which concern methods of the designing, manufacturing and operating processes of such technical objects 
and devices as : ships, port equipment, ocean engineering units, underwater vehicles and equipment as well 


as harbour facilities, with accounting for marine environment protection. 


The Editors of POLISH MARITIME RESEARCH make also efforts to present problems dealing with 
education of engineers and scientific and teaching personnel. As a rule, the basic papers are supplemented 
by information on conferences , important scientific events as well as cooperation in carrying out interna- 


tional scientific research projects. 
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MARINE ENGINEERING 


Semi-Markovian models of the process 
of technical state changes of technical objects 


Jerzy Girtler 
Gdansk University of Technology 


ABSTRACT 


The most important problem of operation of all technical objects (devices) is the problem 
of rational (optimum) decision control of the process of technical state changes of the 
objects. Such control can be realized when a model of the process of technical state chan- 
ges is applied. In view of this, a formal description of the operational process of technical 
objects was presented in this paper, as well as it was justified that the model of the process 
may be a semi-Markovian process of finite set of states. Examples of 4-and 3-state semi- 
-Markovian process were given as a model of the above mentioned process of state chan- 


ges of any device. Difficulties associated with forming such model, which arise from application of the 

theory of semi-Markovian processes, were indicated. Therefore the proposed model is the semi-Markovian 

process whose values are technical states occurring in the phase of operation of many technical objects of 

practical importance (e.g. combustion piston engines, turbines, displacement and rotary compressors, 
impeller pumps etc). 


Key words: model, technical object, semi-Markovian process, technical state 


INTRODUCTION 


The process of technical state changes of any technical ob- 
ject is the process of occurence of successive technical states 
casually linked during time t, which appear one by one in such 
a way that the successive state occurs after a determined dura- 
tion time interval of the directly preceding state. The duration 
time intervals of particular states are random variables of con- 
tinuous realizations and finite expected values. 

The process of technical state changes of every technical 
object belongs to the most important processes occurring in 
the phase of its operation. Course of the process should be ra- 
tional, i.e. that whose realization results from an assumed opti- 
mization criterion. Such criterion may be e.g. expected value 
of operational cost of a given technical object or instantaneous 
availablity factor. The last criterion is important when in any 
instant of operation of particular technical objects, t, a task can 
be assigned to them for realization of which their reliable ope- 
ration is necessary [7 , 8]. The control which makes the opti- 
mum course ofthe process of occurrence of successive techni- 
cal states possible, can be realized when the model of the pro- 
cess, which allows for application one ofthe decision theories, 
is elaborated. In the case of such technical objects as combu- 
stion piston engines and turbines, positive-displacement and 
axial-flow compressors two theories are of importance : the 
statistical decision theory and the theory of controlled semi- 
-Markovian processes [8, 10, 11, 12, 13]. Recently the theory 
of controlled semi-Markovian processes is more and more often 
applied with success to solving different problems of durabili- 


ty, reliability and decision-based operational control of various 
technical objects. The theory may also find use in solving si- 
milar problems associated with operation of many technical 
objects including those dealing with control of the process of 
technical state changes of the objects. For this reason the semi- 
-Markovian model of the process has been proposed and its 
practical applicability justified. 


PREMISES FOR ELABORATION OF SEMI- 
-MARKOVIAN MODEL OF TECHNICAL 
STATE CHANGES OF OBJECTS 


From the definition of semi-Markovian process it results 
[12, 17, 20, 21] that such process is stochastic and of a discrete 
set of states, and its realizations are functions of constant inter- 
vals (i.e. those having uniform values within operation time 
intervals being random variables), and right-hand continuous 
ones. From the definition it also results that the process is de- 
termined when its initial distribution P,= P(Y(0) = s; jas well 
as its matrix function Q(t) = [Qi] whose elements are probabi- 
lities of transfer of the process from the state s; to the state sj, 
during the time not greater than t (1 ] ; ij = 1, 2,..., K), are 
known. The elements Qij[(i £j ; ij = 1, 2, ... , k)] of the matrix 
Q(t) are non-decreasing functions of the variable ¢ [12 , 17]. 
Non-zero elements of the matrix can be interpreted as follows : 


Qi(t)-P (W(tn-1) = sj , O1 7 Tn < t| W(15) = si} ; 
5(,8,6 9. 5 1,] 7 2,5,4 $ 1$] (1) 
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The semi-Markovian model of real process of occurrence 
of successive technical states of any technical object can be 
formed only when determination of states of the process is 
possible in such a way as to obtain duration time of the exis- 
ting state in the instant t as well as the state available in the 
instant T+; , stochastically independent on the preceding sta- 
tes and their duration times [6 , 12]. Hence in the modelling 
which is supposed to lead to elaboration of a semi-Markovian 
model of the technical state changing process of given techni- 
cal objects the analysis of changes of technical states occurring 
in real operational conditions of the objects should be accoun- 
ted for. 

In the case of such technical objects as self-ignition engi- 
nes, piston compressors, impeller pumps etc it has been obser- 
ved that prediction of technical state changes of the objects can 
be performed if actual technical state and future conditions of 
their operation are known [5, 6, 7]. This fact which simultane- 
ously means that technical state changes of such objects are 
not strictly dependent on their time of operation, can be highli- 
ghted by means of the following hypothesis H1 : 


the technical state changing process of any technical object 
(understood as a random time function whose values are ran- 
dom variables representing existing technical states), which 
occurs in a rational operation system (i.e. in such operation 
system where operational cost calculations are carried out), 
is the process of asymptotically independent values because 
its arbitrary state considered in any instant t, (n7 0,1 ,..., m; 
To € T1 € ...« Tm) significantly depends on the state directly 
preceding it, but not on the states which have occurred ear- 
lier and their duration time intervals. 


It should be observed that the hypothesis does not contain 
any contradictions which could falsify it in a logical sense be- 
fore testing it. 


The consequences of the hypothesis are as follows [5 , 7] : 


2 the probabilities (pi; ; ij ; 1,j € N) of transfer of the state 
changing process of any technical object from any actual 
state of the object, s;, to any next state s; do not depend on 
the states in which the process has been before 

2 the unconditional duration time intervals of the particular 
states s; of the technical state changing process of a techni- 
cal object are the random stochastically independent varia- 
bles (T; ;1€ N) 

2 the duration time intervals of any possible state s; of the 
technical state changing process of a technical object, pro- 
vided that the next state will be one of the remaining states 
s; of the process, are the random stochastically independent 
variables (T; ;izj;ij€ N). 


The above mentioned consequences reveal the probabili- 
stic law of changing the technical states of any technical ob- 
ject. They are not contradictory to each other, and their logical 
veracity is doubtless. Therefore there is no obstacle to consider 
the consequences as one common consequence K1, in order to 
use it for empirical verification of the formulated hypothesis 
H1. Such verification consists in experimental investigation of 
veracity of the enumerated consequences taken as one com- 
mon consequence K1. 

The hypothesis can be supported by the more detailed hy- 
pothesis H2 dealing with all technical objects fitted with tribo- 
logical units (e.g. self-ignition engines, piston compressors etc). 
The process of changing their technical states is that in which 
the duration time intervals of each of the states are random 
variables. Particular realizations of the random variables de- 
pend on many factors, a.o. on wear quantity ofthe tribological 
units of the objects. For the objects in question (e.g. self-igni- 
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tion engines) it was observed that wear of their sliding tribolo- 
gical units is weakly correlated with time [3, 4, 9, 14, 18, 26]. 
The observation is important because serviceability of such 
objects mainly depend on technical state (i.e. wear) of their 
tribological units. This made it possible to predict technical 
state of such objects by taking into account only their current 
state as well as service conditions without accounting for the 
states occurred before. In order to highlight the fact the follo- 
wing hypothesis H2 can be given : 


the state of any sliding tribological unit of any technical 
object, as well as its duration time significantly depend on 
the preceding state but not on those occurred before and 
on their duration time intervals because its load and both 
wear rate and wear increments induced by the load are the 
processes of asymptotically independent values. 


The statement contained in the hypothesis results 
from two obvious facts : 


* thereisa close relationship between the loading of tribolo- 
gical units of different technical objects (e.g. self-ignition 
engines, piston compressors etc) and their wear [16, 18, 
26, 277] 

* ina longer operation period of any technical object (e.g. 
self-ignition engines, piston compressors etc) there is no 
monotically increasing load changes of their tribological 
units, hence the service loading ofthe units can be assumed 
stationary [2, 4, 15, 22, 23, 24, 25]. 


In order to verify the presented hypothesis H2 and to deter- 
mine if it is true, it is necessary to predict the consequences 
whose occurrence is possible to be empirically checked. 

The consequences K2 which can be concluded from the 
hypothesis (with accounting for loading features ofthe objects 
and their sliding tribological units) were presented in [6]. 

Verification of the presented hypothesis H1 and H2 by 
means of experimental testing the veracity of the consequen- 
ces K1 and K2 can be performed with the use of the same 
reasoning methods which were presented in [3, 4, 7]. 


VARIANTS OF SEMI-MARKOVIAN 
MODEL OF TECHNICAL STATE 
CHANGES OF OBJECTS 


From the presented hypotheses H1 and H2 it results that 
models ofthe technical state changing process (W(t) : t20} of 
many technical objects such as e.g. self-ignition engines can 
be the stochastical processes of discrete set of states and conti- 
nuous time being the duration time of distinguished technical 
states of the objects. The considered models of the technical 
state changing process of any technical object can be mathe- 
matically expressed by the functions mapping the set of in- 
stants, T € R+ , into the set of technical states, S. Hence to 
elaborate such model it is necessary to determine a finite set of 
technical state changes of the objects in question. Assuming 
the serviceability of technical objects as the criterion for di- 
stinguishing the states one can distinguish the set of classes 
(subsets) of technical states (shortly called ,,states"), S. The 
below given set can be deemed a set of the states of practical 
operational importance [10] : 


S Sf = 123.4 


which have the following interpretation : 


(2) 


$ s; —the state of full serviceability, i.e. the technical state of 
any technical object, which makes it possible to use the 
object within the whole range of loads to which it was ad- 
justed in the phases of its designing and manufacturing 


$ s;-the state of partial serviceability, i.e. the technical sta- 
te of any technical object, which makes it possible to fulfil 
all its tasks (as in s;) but at lower values of operational 
indices (e.g. at a lower overall efficiency, hence — in the 
case of self-ignition engine — at a greater specific fuel oil 
consumption) 

$ s,- the state of limited unserviceability, i.e. the technical 
state of any technical object, which makes it possible to ful- 
fil only some of its tasks (e.g. such state which — in the case 
of self-ignition engine — precludes it from working in ac- 
cordance with maximum continuous rating characteristics, 
and in the case of piston compressor — precludes it from 
filling the air receiver up to a required air pressure, etc) 

$ s,- the state of full unserviceability, i.e. the technical state 
of any technical object, which precludes the object from 
fulfilling any task of the set of the tasks to which it was 
adjusted in the phases of its designing and manufacturing 
(for instance such state of engine, which precludes it from 
work in accordance with maximum continuous rating cha- 
racteristics. 


Elements of the set S = (s; ; i= 1, 2, 3, 4} are values of the 
process (W(t) : t2 0} composed of the states successively oc- 
curring one by one, s; € S, and being casually related to each 
other. 

In the case of many technical objects (such as self-ignition 
combustion engines, gas turbines, piston compressors, etc,) the 
distinguishing of the states s; € S (i = 1, 2, 3, 4) is so much 
important because their use is crucial when they are in the state 
sı or in the state s2. However in the second case the objects 
should be used for as-short-as-possible period only after which 
they should be renewed to bring them back to the state s;. The 
states are values of the process W(t) : t 2 0}which is fully de- 
termined if its functional matrix is known [9 , 12] : 


Q(t) = [Qi;(t)] (3) 
as well as its initial distribution is given : 
pi = P{W(0) = sil s €S;1=1,2,3,4 (4 


Depending on an assumed strategy of maintaining the 
technical objects in the technical states making it possible to 
realize their tasks, diferent realization variants of the process 
{W(t) : t= 0} may be taken into account. 

In the first variant a technical state change can occur in 
compliance with the graph of technical state changes presen- 
ted in the Figure. 


Pi2 o T5 


Pay > Ta 23 > T23 


P34 ? T34 


Graph of changes of the technical states 
s; E S(i = I, 2, 3, 4) of the process {W(t): t 20} 


The following initial distribution of the process 
{W(t) : t= 0} can be assumed : 
pı =P{W(0) =s1}51 
pi=P{W(0)=si}=0 fo i=2,3,4 © 
and its functional matrix in the following 
form complying with the graph of Figure : 
0 Qo 0 0 
Q(t) 0 Qí) 0 
Q(t) = (6) 


Qu 0 0 Qy 
Qi, (t) 0 0 0 


The matrix (6) represents changes of the states s; € S(1— 1, 
2, 3, 4) of the process {W(t) : t 2 0}. The probabilities P;(j = 1, 
2, 3, 4) of the event that a given technical object will be in the 
states s; € S(1- 1, 2, 3, 4) determine a chance of fulfilling the 
task by the object. It is obvious that the user of every technical 
object is interested in the object to be in the state s; as long as 
possible. The chance of lasting any technical object in this sta- 
te is determined by the probability determined for an appro- 
priately long time of operation. It means that the limiting di- 
stribution of the process should be determined. 

From the semi-Markovian process theory it results [12, 17, 
21] that the probabilities of changes of states of any technical 
object are determined by the probabilities pi; of the Markov's 
chain {W(t,) : n 7 0, 1, 2,...} introduced into the process {W(t) : 
t2 0j. The probabilities form the following matrix of transfer 


probabilities : - 
P= py ¢ 1 1,2:5:4] (7) 
where : 
pij = P(Wt«) = 5j | W(tn) = si) = lim Q; C) 


The matrix (7) makes determining the limiting distribution 
of the process {W(t) : t = 0 possible. From the matrix (6) it 
results that the matrix (7) has the following form : 


0 1 0 0 
Pa O p 0 
p =| P2 23 (8) 
Pa 0 O P3 
1 0 0 0 


From the theorem given in [12] it results that the limiting 
distribution of the considered process exists : 


P = lim P{W(t) =s;} = 


(9) 
= lim P{W(t) 2s;/ W(0) 2 s;] 
t—o0 j 
which is determined by the expression : 
T:E(T.) 
= J J "e 

Lea , J=1,2;3,4 (10) 

Yn EC) 

k-l 


and the limiting distribution Tj g =1, 2, 3, 4) 
of the introduced Markov’s chain {W(T,) : n= 0, 1, 2,...} 
fulfils the following equations : 


0 1 0 0 
Pa 0 px3 0 

153555754] 0 0 = usas nal 
31 P34 

1 0 0 0 (11) 


POLISH MARITIME RESEARCH, No 4/2004 5 


MARINE ENGINEERING 


MARINE ENGINEERING 


T1 + T2 + T3 + T4 = 1 


Making use of the relationships (10) + (12) 
one obtains the following formulae : 


P;-E(T)M' ; P5- E(TM 
Í -1 
P34—po3E(T3)M  ; P4= p23p34E(T4)M 


and : 


M = E(T1) + E(T2) + p23E(T3) + p23p34E(T4) 


where : 


(12) 


(13) 


E(Tj) - expected value of duration time of the state 
ges (= 1,2,3,4) 

pj  - probability of transfer of the process 
{W(t) : t2 0} from the state s; to the state s; 
($,S9€9 ; 1]71,2, 3,4112]). 


The particular probabilities P;(j = 1, 2, 3, 4) given 
by the formulae (13) have the following interpretation : 


P,; = lim P(W(t) sj) ; Po= lim P(W(t) = s2} 
t—oo t—oo 

P3= lim P{W(t)=s3} ; P47 lim P(W(t) = s4} 
t—oo t—oo 


In the presented variant, such situations are accounted for 
in which the user is able to risk to fulfil his task in the state s? 
of the technical object and even to risk to fulfil some tasks in 
the state s3 of the object. 

The second variant deals with the case when, within an as- 
sumed operational strategy of technical objects, the distingui- 
shing between the states s, and s? is of no importance. Then it 
is possible to consider the simpler process W(t) : t = 0} of 
technical state changes of the objects, namely the model ha- 
ving the set of states : 

S = (s1, $2, $3} (14) 
as well as their following interpretation [6, 7, 8] : 


* the state of full serviceability, sı, which makes it possible 
to use the object in any conditions and in any range of lo- 
ads, to which it was adjusted in the phases of its designing 
and manufacturing 

* the state of partial serviceability, s2, which makes it possi- 
ble to use the object in limited conditions and in a range of 
loads lower than those to which it was adjusted in the pha- 
ses of its designing and manufacturing 

* the state of unserviceability, sz, which does not make it 
possible to use the object in accordance with the purpose it 
was intended for (e.g due to its failure, performing mainte- 
nance operations on its subassemblies etc). 


A graph of technical state changes of such process W(t) : 
t = 0} of objects of the kind as well as its example course were 
presented in [6]. Hence the process is the three-state one of 
continuous realizations (i.e. time-continuous). It may be assu- 
med that if anyone of the states s2 or s3 does not occur then 
every technical object remains in the state sj. And, the set of 
the technical states S = (51, $2, s3} can be considered as the set 
of values of the stochastic process (W(t) : t2 0} of realizations 
constant within intervals, and being right-hand continuous. 

The initial distribution of the considered process 
having transfer graph [6 , 8] is given by the formula : 


1 for 1=1 
P; = PEW) = 8} = 0 for i= 2.3 


and its functional matrix, if the function Q32(t) 
is different from zero (Q32(t) # 0), is the following : 


(15) 
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0 Qin(t) Q(t) 
Q(t) = Q2 (t) 0 Q>3 (t) 
Q3)(t) Q(t) 0 


In operational practice of technical objects any partial re- 
newal strategy of their technical states should not be applied. 
Use of such strategy may lead to occurrence of the unservice- 
ability state s3 and, as a result, — e.g. to occurrence of a high 
economical loss. It means that renewals can be applied in spe- 
cial cases only, e.g. those resulting from impossibility to dis- 
continue realization of the undertaken task and to complete the 
technical state renewal of the object. Therefore operation pro- 
cess of every technical object may be rational when the func- 
tion Q3>(t) = 0. In this case the matrix (16) obtains the form : 


0 Q(t) Qj(t) 
Q(t) = Q;(t) 0 Q33 (t) 
Q;,(t) 0 0 


And, like in the case of the earlier considered processes, 
for the presented process {W(t) : t 2 0) having the functional 
matrix given by (17), the following limiting distribution can 


be determined : 
P, = ECE) pps p;i;E(T5) 
H H (18) 
_ (1 -PPa MET ) 
H 


and: 


H = E(Tj) + pi5E(T7) + (1 7 pi; p2))ECIT3) 


where : 


(16) 


(17) 


P, 


P,, P5, P - probabilities of the event that a given techni- 

cal object is in the state : s1, S2, $3, respectively 

Pi - probability of transfer of the process {W(t) : 
t 2 0j from the state s; to the state s; 

- expected value of duration time of the state s;. 
The particular probabilities Pj = 1, 2, 3, 4) 
given by (18) can be interpreted as follows : 


P, = lim P(W(t) sj) ; Po= lim P(W(t) = s2} 
to too 
P3= lim P(W(t) = s3} 
t—oo 


E(Tj) 


The presented probabilities which determine possible oc- 
currence of particular states of any technical object, are essen- 
tial in making operational decisions [1, 7, 8, 10, 13]. 


SUMMARY 


» Inthe presented considerations it is shown that the process 
of technical state changes of different technical objects (self- 
-ignition engines, gas turbines, piston compressors, impel- 
ler pumps etc) is a process continuous over time and states. 
Because technical state of every technical object is subject 
to continuous changes hence it is possible to consider the 
sets composed of infinite number of technical states. 

> Identification of all technical states of any technical object 
is neither possible nor purposeful for both technical and 
economical reasons. Therefore a need arises to split the set 
of technical states into a finite number of classes (subsets) 
of technical states. 


> Assuming the serviceability of technical objects to be 


a splitting criterion one can distinguish (in the simplest case) 
the following classes of of their technical states : the full 
serviceability state s;, the partial serviceability state s2, and 
the unserviceability state s3. 
The set of the states S = (s;, S2, s] may be considered as 
the set of values of the simplest stochastic process {W(t) : 
t2 0) whose realizations are constant within intervals and 
right-hand continuous. The process, a model of real pro- 
cess of technical state changes of objects, is mathematical- 
ly described by a function mapping the set ofthe instants T 
(T € R,) into the set of technical states S. Elaboration of 
such a model adequate to the process of technical state chan- 
ges of the objects is indispensable for rational control of 
the process. 
The model presented in the form of the stochastic process 
(W(t) : t 2 0} is the simplest and fulfilling the two condi- 
tions : 
e it works like the original, i.e. it realizes analogical func- 
tions 
e it makes it possible to reveal — on the basis of analysis of 
its structure and functioning mode — the new, hidden fea- 
tures of the processes of technical state changes of the 
objects in question, which are represented by this model, 
namely : the reliability indices expressed by the probabi- 
lities, P;, of lasting the process in distinguished states. 
From the presented hypotheses it results that the process of 
technical state changes of the objects in question can be 
investigated by means of the models formed as semi-Mar- 
kovian processes. Therefore the process (W(t) : t2 0} can 
be considered as a semi-Markovian process of real proces- 
ses of technical state changes of many technical objects. 
The hypotheses should be verified for each kind of techni- 
cal objects, e.g. self-ignition engines, piston compressors, 
impeller pumps etc, as it may happen that for a given kind 
of technical objects operating in assumed service condi- 
tions it would not be possible to use the models of technical 
state changes in the form of semi-Markovian processes. 
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Miscellanea 


Flutter of Turbine Rotor Blades 
in Inviscid Flow 


In the year 2004 Polish Naval University in Gdynia edi- 
ted (in English) the book on : "Flutter of Turbine Rotor Bla- 
des in Inviscid Flow" 


written by Prof. Romuald Rzadkowski 


The motivation to write this monography was to present 
theoretical models, numerical codes and results of the aero- 
elastic calculations of rotor blades which have been done 
since 1997 in the Institute of Fluid-Flow Machinery, Polish 
Academy of Sciences in Gdansk. 

These results were obtained in collaboration with Prof. 
V. Gnesin from Institute for Machinery Problems, Ukra- 
inian National Academy of Sciences, Kharkov branch, and 
dr V. Tsimbalyuk from the Institute for Strength Problems, 
Ukrainian National Academy of Sciences, Kiev. 

Turbine blading dynamics is an interdisciplinary doma- 
in of knowledge, which inter-connects mechanics of mate- 
rials, vibration theory, aerodynamics, thermal exchange, pro- 
blems of corrosion and erosion effects of surrounding me- 
dium, and manufacturing technology. The most important 
information required to perform the dynamic calculations is 
the blade unsteady force which arises out of the flow path 
interaction in the stage and the flutter parameters. 


The presented book is limited to the flutter 
of rotor blades in inviscid flow. 


It includes physical and mathematical modelling of self- 
-excited vibration in 2D and 3D flows. For building mathe- 
matical models of blade and bladed disc motion the general 
theory of continuous media and methods of analytical me- 
chanics were employed. 


The entire scope of the book 
is contained in the seven chapters. 


O Chapter 1 is introductory 


O In Chapter 2 the flutter phenomena and aeroelasticity 
methods are described. 


O Chapter 3 presents the 1D and 3D structure models of 
a blade and a bladed disc. 


O In Chapter 4 the 2D inviscid flutter model is given. The 
2D trans-sonic flow of an ideal gas through a multipas- 
sage blade row is considered. In a general case the flow 
is assumed to be an aperiodic function from blade to bla- 
de (in the pitchwise direction), so the calculated domain 
includes all blades of the whole assembly. The aerody- 
namic model fully accounts for the blade thickness, cam- 
ber and the angle-of-attack effects. The unsteady flow of 
the ideal gas is described by the 2D Euler equations. 
Godunov-Kolgan method is used to descritize the 2D 
Euler equations. The numerical test calculations were per- 
formed to compare the theoretical results obtained by 
using mathematical model developed here, with experi- 
mental data. 


O In Chapter 5 the model of 3D self-excited vibration of 
bladed disc in the compressible flow is shown. The 3D 
non-linear time-marching method for aeroelastic beha- 
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viour of oscillating turbine blade row was presented. The 
numerical calculations were performed to compare the 
theoretical results of 3D inviscid flutter code with expe- 
riments. The aeroelastic behaviour of a 3D oscillating 
row of last stage turbine blades of L=0.765 [m] in length 
is shown for a harmonic motion and coupled fluid-struc- 
ture interaction. 


O In Chapter 6 the experimental stand to simultaneously 
measure unsteady aerodynamic force and moment with 
arbitrary combinations of blade motions in the subsonic 
flow, is presented. The numerical calculations were per- 
formed to compare the theoretical results with experi- 
ments for the harmonic motion. 


O Chapter 7 presents the general conclusions and remarks 
and the directions of the future work. 


The book is addressed to the scientific workers and stu- 
dents of the university faculties engaged in the field of aero- 
elasticity applied to steam and gas turbines. 


The presented book was preceded by another interesting 
book of the same author, titled: Dynamics of Rotor Steam 
Turbine Blading, Part Two, Bladed Discs edited by Ma- 
szyny Przeplywowe, Wroclaw, Ossolineum, 1998. 


Prof. Romuald Rzadkowski graduated from Shipbuilding 
Institute, Gdansk Technical University, and Faculty of Ma- 
thematics, Gdansk University. Presently he works in the In- 
stitute of Fluid-Flow Machinery, Polish Academy of Scien- 
ces, and Polish Naval Academy in Gdynia. He obtained his 
PhD degree in 1988 and DSc in 1998 at the Institute of Fluid- 
-Flow Machinery, Polish Academy of Sciences. 

Prof. Romuald Rzadkowski is the author and co-author 
of about 123 scientific papers. The area of his scientific inte- 
rest covers free and forced vibration of turbine blades, bla- 
ded discs, shaft with bladed discs, life estimation, flutter and 
unsteady forces in the turbine stages. 
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Viscoelastic lubrication of spherical slide 
bearing in impulsive unsteady motion 


Krzysztof Ch. Wierzcholski 
Gdansk University of Technology 
Gdynia Maritime University 


ABSTRACT 


This paper presents numerical and semi-analytical solutions of oil velocity components 
and pressure distributions in spherical unsymmetrical gap of slide bearing. A hydrodyna- 
mic unsteady lubrication during oil flow with viscoelastic properties is here considered. In 
the case of various driving systems on ships the bearings with spherical journals and 
spherical sleeves or slide bearings with spherical bits operate often under impulsive unsteady 
motions. Many impurities appearing in service leads to viscoelastic properties of the oil. 
During service of transport machines it is necessary to adjust the shaft location respective 


to the sleeve in order to make optimizing the convergent lubricating film possible. Such conditions are 
effectively satisfied in bearings with spherical journals. The presented numerical calculations were perfor- 
med by means of the Mathcad 2000 Professional Program and the method of finite differences. This method 
satisfies stability conditions of numerical solutions of capacity forces occurring in spherical bearings. 


Key words : driving systems on ships, spherical slide bearing, unsteady impulsive viscoelastic lubrication 


INTRODUCTION 


Lubrication of spherical bearing under periodic motion has 
been considered in many papers till now [3, 4, 8, 9, 10, 11, 12]. 
This paper considers pressure distribution during hydrodyna- 
mic viscoelastic lubrication of spherical bearings and bearings 
with spherical bit, at impulsive unsteady motion. These pro- 
blems have been not elaborated hitherto. 

Bearing systems are commonly used in diesel engines in- 
stalled in land transport machines and ships. The oil in bearing 
gaps in such bearing systems is contaminated mainly with dust, 
soot, smoke black as well as many inhibitors improving the oil 
properties. Transport machines usually work under unsteady 
impulsive vibrations. Thus the lubricating oil has often the non- 
-Newtonian properties. Therefore in this paper viscoelastic time- 
-dependent properties of oil are taken into account. Designing 
the bearings without accounting for the viscoelastic oil proper- 
ties brings about to occurrence of the seizing of the bearing in 
kinematics pairs. 

The seizing of bearings can be prevented by proper reco- 
gnition of bearing operational parameters for real oils. This is 
very important because the seizing of ship diesel engine bea- 
ring system may lead to the catastrophe [1, 2]. Therefore deter- 
mination of real bearing capacity at unsteady viscoelastic lu- 
brication has important sense. 

The bearings presented in this paper have spherical jour- 
nals. The spherical journals can be turned with the shaft or the 
form of an individual ball can be used (see Fig.1). Such ball is 
installed in bored end of the shaft. 


Fig.1. Spherical bearings: a) spherical journal together with shaft and 
spherical sleeve, b) individual spherical ball with conical sleeve, c) spheri- 
cal journal together with shaft and conical sleeve 


The sleeve has spherical or conical shape. The spherical 
sleeve is more effective than conical one, because it ensures 
small values of slide thrust and wear. 

Usually the spherical bearing is adjustable. It makes it po- 
ssible to set up the shaft with respect to the sleeve and to con- 
trol the convergent lubricating film, and it is capable of trans- 
ferring both axial and transverse loading, Fig.2. 


Fig. 2. Loading of the spherical journal: a) transverse loading, 
b) axial and transverse loading 
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GOVERNING EQUATIONS AND BOUNDARY LAYER SIMPLIFICATIONS 


Lubrication of journal and sleeve in slide spherical bearing is described by oil flow. The momentum conservation equations 
and continuity equation describe oil flow. Moreover the second order approximation of the general constitutive equation given 


by Rivlin and Ericksen, can be considered. The equations can be expressed in the following form [6 , 7] : 


Div S- pdv/dt , divv=0 , S=-pl+noA1+ o(A4). + BA? 


where : 
S - Stress tensor p -oil density 
Div S - stress tensor divergence t -time [s] 
v - velocity vector [m/s] p -pressure 
divv - velocity vector divergence I] - unit tensor 


A; and Ay - two Rivlin-Ericksen strain tensors of three material constants No, a, p, 
where : 
No - dynamic viscosity a , D - pseudo-viscosity constans of oil. 


The tensors Ay, and A»? are given by the symmetric matrices defined by: 


Ay = L sr E 
A» = grad a + (grad a)! +2L'L 

we 

a=LVT =Z 

ot 

where : 
L -tensor of oil velocity vector gradient [s] a - acceleration vector [n/s?] 

Ll -tensor with matrix transpose [s ] grada -acceleration vector gradient 


(1) 


Q) 


It is assumed that the product of Deborah and Strouhal numbers, i.e. DeStr, and the product of Reynolds number, dimension- 
less radial clearance, and Strouhal number, i.e. RewStr , are of values of the same order. Moreover DeStr >> De = o/n, . 


where: w -relative radial clearance € - angular velocity of spherical bearing journal. 
The following is additionally assumed : 


the rotational motion of spherical journal with peri- 
pheral tangential velocity U = œR 

unsymmetrical unsteady oil flow in the gap 
viscoelastic and unsteady properties of oil 

the oil density p of constant value 

the characteristic value of the bearing gap height, € 
no slip at the bearing surfaces 

R - radius of spherical journal. 


VVVVVV WV 


By neglecting the terms of the radial clearance y = €/R = 10° in the governing equations expressed in the spherical 
coordinates @, r, 0, and by taking into account the above mentioned assumptions the following is obtained : 


3 
Ae LV MSIE Yo 
2 
pRsin2.00 P orl Or J p ator 
R 
pe 
Or 


an (Aii 0 vs 


Ov 
— + Rsi >| + a Rv, LE -0 
Co) R/ or ð R 


where: 05 px2mc; , O«c4«1 , bn=nR/8 <8 <nR2=b, , 0<r<h , h- gap height. 


(3) 


(4) 


(5) 


(6) 


Symbols Vo, Vr, Vg denote oil velocity components in the circumferential, gap-height and meridianal directions of the sphe- 
rical journal, respectively. The terms multiplied by the coefficient B in the right hand sides of (3) , (5) denote influence of time- 
variable viscoelastic oil properties on the bearing operational parameters. The terms in the left hand sides of (3) , (5) describe 
influence of accelerations which occur during the impulsive motion, on the bearing lubrication. The relationships between 


dimensional and dimensionless quantities are assumed in the following form : 
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r-gr , SROI , t=toty , h=ehy , vg=UVgi , Yr= Uyvg 
2 
Vo =Uve1 , P=PoP1 . pos UNoR/E 


and Reynolds number, modified Reynolds number and Strouhal - Deborah number are as follows : 


Re =pUe/n, , Rey = POE? Mo , Str = R/(Uto) , Des BU/(MoR) 
hence : 


DeStr = B/Moto) = Des , ReyStr = pe / Moto) = Res 


(7) 


(8) 
(8a) 


For the oil containing inhibitors the constant B/t, is always 0 < B/to < No and of values usually in the range from 0.0001 to 
0.1000 Pas”. The dimensionless symbols have lower index “1”. Hence the equations (1) + (6) obtain the following dimensionless 


form : 
3 
Res ONoi "PEE AA pi, O0. [27H + Des LZ 
0 = Op, 
Or, 
Ras ÓVgi E Op; ð | OVg " ET 
at, 09, On Or at, ar, 


where:0€ox2mc, , OSc) <1 , m/8x804 m2 , OSr Shy 
THE METHOD OF INTEGRATION 


For lubrication at impulsive motion a new dimensionless variable was introduced [4] : 


I IR D 
y*nN, s E , t,>0, paa] 


t 
1 1 
and the solutions in the form of the following convergent series, were assumed : 


2 
Des Des 
Vol = V o0x (X. 9. 91) + FED (x. 9. 91) t (Pe) V 92 (x. Q9. 91) Teen 


1 1 


2 
Des Des 
Vai 7 Vgoz 09. 94) + PEL» (59.91) + e V 2x (Xo 9. 91) t ......... 


1 1 


Des 


Des 
Vit = Vros p P, 91) + t E 


2 
d 1 | Vix 059, 94) + ein 


1 1 


Des 


Des 
Pı -pu (dot) + pesa] t | Pp (0, 9,, t, ) + — 
1 


1 


where: t; >0 , 0<Des<<1 , (Des/t})<1 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


The first terms of the series (14), (17) describe oil flow parameters in impulsive unsteady motion, with neglecting the visco- 
elastic properties. The second, third, etc terms in the series (14), (17) describe the corrections of oil flow parameters, caused by 
the time-changeable viscoelastic oil properties. In (9) - (11) the derivatives with respect to the variables t;, rj, can be replaced by 


the derivatives with respect to the variable x only, by using the following relationships : 


8 (0)_a( 8 oy) ax Res e 
Or? Or, Or, ) Oy Oy Or, J Orn, 4t, Oy? 
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e? ô (Res &? Res 0? Res O( 0? \ ay Res(0) y% æ 
jw e aaa a a oo le a a a (19) 
Ot Ore ot, | 4t, Oy At; Oy 4t, OX Oy jat 4ti (Oy 2a 

Next, the series (14)-(17) were put into the changed set of the equations (9)-(12) where the variables t;,rj were replaced by the 


variable x. And, the terms multiplied by the parameter (Des/t1) ` of the same power values k , for k = 0,1,2, ..., were respectively 
equalled to each other.Thus the following sequence of the sets of ordinary second-order differential equations, was obtained [5] : 


(Vios r + 2X(Viox j" p Im (20) 
Q) D _ 1 Pu a, 1 (2) "T 
(Vis) tme) das) — TV JE (s^ (21) 
N? 0a; 2 
(2) () _ 1 Opp Q1 3) 
(Vios *2X(Viz) *8(vi5x)-7 NI zm *t2(Vug)"- 2 XC ux) (22) 


i i 
and so on, 
where:i79,0 , Ap=Q , 095204 
The upper indices: (1), (2), (3), ... denote : the first, second, third, etc. derivative with respect to the variable x, and : 


(No) =N’sin(0;) , Ng=N (23) 
GENERAL SOLUTIONS AND VALIDITY OF BOUNDARY CONDITIONS 
The general solutions of the equations (20) for : i= @ , 0; have the form : 
Viox(X) = Ciivo100 + Ci2voalX) + vi 03(X) (24) 
where : Cj; , Ci? - integration constants. 


The following particular solutions of homogeneous and non - homogeneous differential equations were obtained : 


X 2 
vy GO =e" dy, , vS 0071 (25) 
0 
1 Op X, 2 X, 2 
Vig; (X) 7 - a p” Voici dx, — Vor) je" dy, (26) 
i i13 0 


where : 0 € X1 € X =r1N. 


For ti > 0 ,N — œ , thus X, > os. For t; — 99, N — 0 hence for rı > 0 will be x, > 0. 
For t; > 0 and rı = 0 also x = 0. The following limits are true : 
voi(x) = n9? for:y 05, tj 0, No 
Vo1(X) = 0 for: y% > 0 ,r170 , 0<ti<t<o , N>O 
Vo1(X) = 0 fon:y0 , n»0, to , NO 
Vios(X) = 0 fon:y 0 , m=0 , 0<tı<tı «o , N»O where:i=9,9 (7) 


2 
r Opi 
M =- ——-  fo:*y0,r7»50,t,59 ,N—0 
s Qo 2sin8, Op X i i 
2 
vas) =- Po fo:Yy0,rT,»0,t, 59 , N>0 
2 09, 


The spherical journal moves only in the circumferential direction @. Hence the oil velocity components on the journal surface 
in this direction are equal to the peripheral velocity of the spherical journal surface. The oil velocity component on the spherical 
journal surface in the meridianal direction 9 equals zero because the spherical journal is motionless in 0- direction. The oil flow 
around the journal is assumed viscous. Hence on the journal surface the oil velocity component in the gap height direction equals 
zero. Therefore the following boundary conditions are valid: 


Veox(y 7 0) = sin81 , vsox(4 70)7 0 , vix 70) =0 a 
fo:rj-0 €» y -OandO «tj «t)5«o , N>0 
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The spherical sleeve surface is motionless in both circumferenial and meridianal directions. But the spherical sleeve performs 
any impulsive displacements in the gap height direction. Hence the gap height changes along with time. Thus the oil velocity 
components on the sleeve surface are equal to zero in both circumferential and meridianal directions. The oil velocity component 
in the gap height direction r is equal to the first derivative of the gap height with respect to time. Hence the following boundary 
conditions are valid : 


Vooz(, =M)=0 , vsox(x, 7 M) * 0 , vosí% = M) = Strohj/ot; 

(29) 
for : Ty >h ey% —>Nh =M , 0<t<t «o , N>O 
where: h=eh, - gap height , hı - dimensionless gap height , Str = l/ot; 
Imposing the conditions (28) , (29) on the solution (24) one obtains : 
Co1Vo1(X% -0)- Co2 T Voo3 X, -0)- sin9,  for:ri-0 
Co1Vo1(X, = M) + Co + voos(, = M) = 0 for : r1 = hy 
Cg1Vvo1(% = 0) + Cgo + vgos(y, =0)=0 for: r= 0 
Ca1v01(% = M) + Cy + Vgos(y, = M) = 0 for: r1 = hy 


(30) 


By taking into account the limits (27) the following solutions of the set of the equations (30), are obtained : 


sin 9, + v M V M i 
"E 1 mt ) om —_ 303( ) , Co =sin 9, ; Ce E (31) 
Vo, (M) Vol 
Now, into the right hand side of (21) the solution (24), (25), (26), (31) is inserted. 
Thus the general solution of (21) obtains the following form : 


Vas) = Cay O + CaviOD + Vis) for: 1 =O, 8 (32) 
where : Ci; , Ci4 - integration constants. 


The particular solutions are as follows : 


y 2x1 E 
vu) =xe* , vo) xe* ae “dy (23) 
ô X1 


Viso Ci) = vi) Ci Xx * 2) zm A je Le ri Vi03 Qu] E Vio (%ı)dxı + 
i l (34) 
+ Viz col t F A je y FI (x)0]*— N? c: 


fot:i29,80 , 0<d<4%, <x 


Cii 0G + > te 


The solutions (32) represent the corrections of the oil velocity components due to the viscoelastic oil properties. 
By virtue of the solutions (33) and (34), for: y 0 , rj 90 , N>O, it follows: 


lim Vi(X) = lim xe * [7 eMdy, = (35) 
y0,N»50 y0,N»0 PA 
The following limits are true : 
vii(x) - 0 fot:y0 , m=0 , O«ti«t «o , N>0 
Vix) = -1 fo:y 0 , n=0 , O«tj«t «o , N>0 (36) 
Vis(X)70  fo:y 0 , m=0 , O<ti<t<o , N>O whee:i=ọ,9 


The corrections of oil velocity components can not violate the boundary conditions (28),(29) which are assumed on the 
journal and sleeve surfaces in the circumferential and meridianal directions. Therefore, the following boundary conditions were 
applied to the corrections of oil velocity components : 


Von ( —0)— 0: , was(y-0)70 Yoe:r—0c»y-0 y Uisp ,N-—U 


(37) 
vaz% 7 M)7 0 , Vsis(y=M)=0 for:r, hie y —5NhieM , 0<ti<xtp<o , N>0 
Imposing conditions (37) on the general solution (32) one gets : 
Cavu = 0) + Cava = 0) + vois( 70) 70 — for:r;-0 ES 


Cos iX x M) + Co4V21(%, = M) + Vo13(X = M) =() for: T1— hi 
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Co3vi1(%, = 0) + Cava% = 0) + vgs% = 0) = 0 for: r= 0 


(39) 
CasVii (X = M) + Ca4Voi(XX = M) + Vois =M)=0 for: r= hi 
By taking into account the limits (36), the following solutions of the set of the equations (38),(39) were obtained : 
MEUS LUE Cií 20  fo:i-209,8 (40) 
vix = h,N) 
where: 0S xj £hiN , N= ee | OK<t<0 , O<t Shy , bai $B, <be , O<O<2mc, , Osce 
1 


NEWTONIAN UNSTEADY LUBRICATION 


By neglecting the viscoelastic properties of oil, by virtue of solutions (24) and constants (31), the particular velocity 
components of oil in @ - and 9 - directions for non steady flow obtained the following dimensionless form : 


Voos (9.15,9,.t,) = sind, — sin9, - vn PP Fy (y = Nh, J . erf(r, N) " 


vn ð 
—3— — "PI y (y= Nr) 
2N? sing, Op 
vn Op erf(r,N) vn Op 
V n,9,,t,) 2 — 2 lY (y = Nh, | ——ÀÀ—- ——2 Y(y =N (42) 
sox (1.9.1) IN? 09, | (x i) erf(h,N) 2N? 09, (x ) 
where : 
X 2 ES 
Y(y)s fe" erfy dy, - erf(h;N) fe? dy, (43) 
0 0 
1 [Res 2 X .2 
1l NU UE (44) 
7 b (ua) Tee X2 
for : 


O0<t<œ , OSrjEh, , bm SA Sbs , O<QM<2mc] , OScy<—~ , O<X2<SX1 SX=r1N <h N=M , hı =h (9,®t1) 


The oil velocity components (41),(42) were put into the continuity equation (12) and both sides of this equation were integra- 
ted with respect to the variable rı. The oil velocity component v;ox in the gap height direction equals zero on the spherical journal 
surface. Therefore by imposing the boundary condition v,gs = 0 for rı = 0 , the oil velocity component in the gap height direction 
obtained the following form : 


Ne" Jah; Vaf 1 2h pu, hi io x feet | OY a, ș 
erf(h,N)| 0p 2 |sin?9, dp ao * 69, 08, perf(h,N) ^ 
D o e? ð 1 " erf(rN i 

E > Pu Ho , "Pio cot g, — Y( S eese) dr, - [Y(x 2r, N)dr, 
2 isin^89, OQ O9; 08, N oerf(h,N) 0 


where :0Stj<0 , O€rosr sh; , by SO Sb, , O<QM<2mcy , 0xcqj «os , OS yoSy1S X SrINShjIN2SM 


Vos (9.15, 9.1) 2 — 
(45) 


The oil velocity component v;ox in the gap height direction does not equal zero on the sleeve surface. Therefore by integrating 
the continuity equation (12) and imposing the boundary condition (29) for r;-h; on the velocity component in the gap height 
direction, the following equation was obtained : 


h h ah 
ee dr, + = Í sin 91V goz dr, = “Str sin, (46) 
P o 1 0 1 


If the expressions (41)+(42) are put into (46) the following modified Reynolds equation is yielded : 


hy 
G 4 8 l erf(r,N)dr, 
2N? sin9, 0o ||  erf(h,N) 


_ M _ OP 9 
Y(x = Nh, ) J Y(x = Nr, )dr, a zd 
0 (p 
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hy 
i us l erf (r, N )dr, h 


Op 
—_ — 5] 1 Nh Nr, Jd 
+n? 29,11 efN) Y(z=Nh,)- [Y (x — Nr, dr, 25. 


9. sin 4s (47) 


MARINE ENGINEERING 


h; 
- - (sin 9, )— i [ jou dr, agen 
OP | o erf(h,N) Ot, 


where : 
O€rtzrj€h,, OxXqo«2mce,, OScy <1 , 0x0, «m2 , Oti «eo 
O<% <% ShN , 0<N(t1) =0.5(Res/ty)°? < 
The modified Reynolds equation (47) determines an unknown pressure function p19 (@, ®1,t1). If tı tends to infinity, i.e. N 


tends to zero, then the equation (47) tends to the classical Reynolds equation. To explain this fact the following limits were 
calculated : 


TU 
Be esi Res 
Mn NP (ce hN) 


h,N hN 
i | | exp GC Jerf(x)dy — erf (hN) l exp oa = 


1 |hN T " h,N " h,N : 
7lim—51! E )Jexp xi s lx- f exp (-x ve] f exp(x Lr] - (48) 
N50 N 0 0 0 0 
Nh, , 
h 
H | expo dx H ph? exp(h?N?) 7 h? 


- lim ——— ———— - - — [im ———————À—————----— 
N50 2Nexp(h?N?) 2 NoO0exp(h?N?) - 2h? N?exp (h2N?) 2 


and, analogously : 


2 
Vn r 
lim —— Y(x = Nr, )= - — (49) 
cae IN? (x 1 ) 2 
as well as : 
erf(nN) m (50) 
N30 erf(h,N) ^h 
Thus the equation (46) for N > 0 tends to the following form : 
1 ð hy )^ h ô h? 
- —4||-— far —f a dr, PD re | oa sin 9, (Lar, + 
sin 9, 0) 2 0 h, 0 2 op “39, 2 oh, 
(51) 
TE x ô sik oh 
E iL sin 9,dr, Pig = —(sin 9, )— | o dr, |- Str — (sin 8,) 
$V 3 09, Op|oN h; ot; 
Finally, after calculations, the following form of the classical Reynolds equations 
of flow in the spherical coordinates was obtained : 
1 Ó Oh Oh 
ni Pro |, li Pro si in9, |- 6— sin 9, +12Str—sin9, (52) 
sin 9, E, OQ taS a3, OQ et, 
lor:0sS9p«2mc; , 0<cı <1 , 0801-02 
The time-dependent average gap height with perturbations has the following form : 
= (ho/€)[1 + sı * exp (-tot100)] 
where : (92) 
ho(@,81) = A£,cosqQ sin, + A£» sing sind; — A£3 cos0; - R + (54) 
0.5 


t [(A£1c0sqQ sind; + A£» sing sind; — A£3 cos0,) + (R + Eun) (R + 2D + £u] 
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The coefficient s controls changes of the gap height during the impulsive motion. If s > 0 the gap height increases, ifs < 0, the 
gap height decreases. The symbol c denotes an angular velocity expressed in [s ] , describing impulsive changes of perturba- 
tions in the unsteady oil flow in the bearing gap in its height direction. If Strouhal number tends to zero the equation (52) tends 
to the classical Reynolds equation for stationary flow. 

The centre of the spherical journal was assumed in the point O(0,0,0) and the centre of the spherical sleeve in the point O, 
(x — Ag, y — A£5, z + A€3). The eccentricity was determined by the value D (Fig.3). The lubrication region Q indicated in Fig.3 
was defined as follows : 0 € o < n, nR/8 € o3 = 0 € tR/2. 


; force of impulse 
spherical sleeve 


Fig.3. Schematic diagram of the gap height and eccentricities 
VISCOELASTIC UNSTEADY EFFECTS 


The particular corrections (32) of the oil velocity components in @ — and 94 — directions, caused by the viscoelastic, fluid 
properties and unsteady fluid flow, were multiplied by the factor DeStr/t;. By using the expressions (32), (33), (34), (40) and 
boundary conditions (37), the corrections of the oil velocity components (32) obtained the following form : 


Des 4B Nre Op 
Voix (0.94.1. t )2— p Seia L'i 
op 


Yi (x) dy + 
us pe? sin 9, Ix (Ddy 


nN 


yf 2 2) u , Pro |" h y d 
m" h2- r? JY, (; 2 h,N) |+ o TY, GO dx - Y (x = ND fe Y> (%)d% + 
nN 
55 
ô Pio d 


Jnerf(h,N) 09 


tN 2 
t Yi(y 2 nN) xe ^ Y (x) ax- h JY 00 Y4 Q0Y Q0 dx + 


hjN 2 
- Yi (x 2 hN) l Y3 (x) ¥(x) dy + Yi 2 nN) Y3(K)Y(X) ul + 


Yx% = hN) T Y3(%) dx - YG=AN) | Y3(%) dx - "TY, GOY; war 


nN 


T BBN re N’sin 3, | 
vnpe’erf (h,N) 


4 N? 
2 v9 Stt )= = No Tale [xYi God + (hi m 2)y coh] 


1 1 LN 
O 
Pio 


1 


_ 2 ô P10 
vnerf(h,N) 09, 


i Y, 60d — Y (x = h,N) i xe * Y () dy + Y, ( = nN) i xe * Y id: (56) 


nN 


ki Yi Y4 Q0 Y (dy, — Y; (x, = h,N) T Y400 YQ dy Y (y =N) j YOY al 


nN 
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with : 


X] . y? 2X 9 2 

= x = = X — S 

You dy, , s) [c 20 s Jerda -1| » Ya(x)=x°(x+2)e* (57) 

6X1 0 
whereas: O Stj <œ , OSm<ry<hy , bm SB) Sb, , 0«qg«2nc; , OScy<~m , OS yi SX SrjINShjINSM 

The corrections of the oil velocity components (55) , (56) were put into the continuity equation (12) and both sides of this 

equation were integrated with respect to the variable rj. From the viscous oil properties it follows that the corrections of the oil 

velocity components in the gap height direction equal zero on the journal surface for rı = 0.Thus the corrections of oil velocity 

components in the gap height direction obtained the form : 


2 ô 
sin 9, Phasen jan |- sin 9, 09, ilins, )vgis (9.9,.5.t, «| (58) 


The corrections of the oil velocity components can not violate the boundary conditions (28) , (29) 
assumed on the journal and sleeve surfaces in the gap height direction. 
Hence the corrections of the oil velocity component in the gap height direction equal zero on the sleeve surface for r = hy. 


Vas (9,9;,1,t;) =- 


By imposing this condition on the solution (58) the following modified Reynolds equation was obtained : 


1 ð LU ‘Tse NS (dr, + Z,(h, JE E zm sin 8 ICE ene) 
0 


sin 9, 09 | 0o 09, | 08 
_ 2Nsin 9," a 2 0 JO : —1 N° 
L fne He rN’ Z3(4) dt; 4 ———— ————— — at Lue E» Z4Q)dr -Zs(h)| a 
~ erf(Nh,) )o V nerf (h;N)sin8, Op | OM Lo (59) 
2 OP io | fie N? 
"EE SM sin9,| [ne * " Z,(r)dr - Z;(h) |) 
* Mnerf(h,N) 09, m ~ — 


hy 272 
D 0 [eno Iz ne s E, s(t, )dr, eZ) | Piana, je Z(t, dr, ez) 
1 


sin, 0p | Op 09, 0 
- where: 
1 ] age 
Zn) = [xY G)dx- n Yi(c- hN) . Z;)- zü-e "hr Yi Gi = Nh) (60) 
rN 
Zin)2 Yi(x = gN) [Ys (Ody - Yi (x = aN) f Y3 (x%)d% - T Yi(X)Y3 (Wy, (61) 
rN 
Z,(n)7 T Y 00 YS CO ¥ dx +Y x = aN) i YiWY3 G0 Y Ody (62) 
rN 2 N? 
—e 
Zs(hi) REC OX (x =h N) ii Y WY, 60 Y (x )dy (63) 
hN 
Z(t = fY ()dx * Y (x = =1N) f AA dy (64) 
r,N 
i= pn" hN , 
Z4(h, = X (x = h,N) ae (x) xe "dx (65) 


and:0<r2<rı <h; , 09 s2neqp , 0Scj«1 , 0S0, 12 , ÜSt( <% 
O<% XLN , OSN(t) = O.5(Res/t)) ^ < 
The modified Reynolds equation (59) determines an unknown function p11(@, ®1,t1) 
of the pressure corrections due to the viscoelastic properties of oil in unsteady conditions. 


NUMERICAL CALCULATIONS 


The dimensionless pressure distribution po and its dimensionless corrections p11, p12 , ... are determined in the lubrication 
region Q by virtue of the modified Reynolds equations (47) , (59) and by taking into account the gap height (53) , (54). On the 
boundary of the region Q the dimensional pressure and its corrections have values of the atmospheric pressure pat. The region Q 
indicated as a section of the bowl of the sphere (Fig.3), is defined by the following inequalities: Q : 0 € p € n, nR/8 € 0 € nR/2. 


Numerical calculations were performed for : 


& the radius of spherical journal R = 0.08 [m] 
© the angular velocity of the perturbations of bearing sleeve œ, = 0.2 [s'!] 
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« the characteristic time t, = 0.001 [s] 
« the characteristic value of the radial clearance V = €/R = 0.001. 


And, the following bearing eccentricities Ae1=20 [um] , Ae2= 2 [um] , A£3= 1[um] , the oil viscosity No = 0.03 [Pas] , the 
pseudoviscosity coefficient B = 0.0006 [Pas?], the oil density p = 950 [kg/m], the rotational velocity of the spherical journal, 
n = 1500 [rev/min], and the average minimum gap height Emin = 4 [um], were assumed. 

The numerical calculations were performed by using the Mathcad 11 Program and the finite difference method. The obtained 
pressure distributions for the dimensionless time values tı = 1, t; = 10, t; = 100, t; = 1000, t; = 10000, tı = œ, and s =+1/4, are 
presented in Fig.4 and 5. 


Rees 
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t=l 
Emin=5.52 um 


Pio may 0.403 
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X o. 26 y 
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-0.13 i 
x p l 26 y 


pio pio 


£i 9.2 0pm 
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-0.26 -0.13 


x Qo — un 


-0.26 -0.13 


X 
Fig.4. The dimensionless hydrodynamic pressure distributions inside the gap of slide spherical bearing, over the region 42: 0 X p <m, nR/8 < 0 < 7R/2, 


at the dimensionless time values : tj = 1, tı = 10, t; = 100, up to the impulse occurrence, for the increasing (decreasing) effects 
of the gap height changes, shown in the right (lefi) hand side column of the diagrams, respectively 
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In order to obtain real values of time it is necessary to multiply the dimensionless values tı by the characteristic time 
to = 0.001 s. For example t, = 1000 denotes the time of 1s after impulse occurrence. In order to obtain realistic dimensional 
pressure values the dimensionless pressure values indicated in Fig.4 and 5 are to be multiplied by the dimensional coefficient 
URNE“. 

The pressure distributions shown on the right-hand sides of Fig.4 and 5 were obtained for the increasing of the gap height, caused 
by the impulse effects. In this case the longer the time up to the impulse, the more gap height decreases and pressure increases. 
The pressure distributions shown on the left-hand side of Fig.4 and 5 were obtained for the decreasing of the gap height, caused 
by the impulse effects. In this case the longer the time up to the impulse, the more gap height increases and pressure decreases. 
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Fig.5 The dimensionless hydrodynamic pressure distributions inside the gap of slide spherical bearing, over the region 2: 0 <o x zt AR/8 $0 < 7R/2, 


at the dimensionless time values: t; = 1000, t; = 10000, t; — ee, up to the impulse occurrence, for the increasing (decreasing) effects 
of the gap height changes, shown in the right (left) hand side column of the diagrams, respectively. 
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If the time interval up to the impulse occurrence is sufficiently large i.e. for tj — œ , then the pressure distributions for the 
increasing (s; > 0) and decreasing (s, < 0) effects of the gap height changes caused by the impulse, tend to the identical pressure 
distributions (Fig.5). Such limit pressure distribution can be also obtained from the classical Reynolds equation (52). 

For the dimensionless time values : tj = 1, t; = 10, t; = 100, t; = 1000, t; = 10 000, t = eo, i.e. for: t = 0.001s, t= 0.010 s, 
t= 0.100 s, t = 1.000 s, t = 10.000 s, t = œ s, after impulse occurrence the maximum pressure distributions for s, < 0 have the 
following dimensional values, respectively : 


111723. , 111525. , 11022. , 0,9932. , 0,6742. 


V V Vy V V y 

For the dimensionless time values : tj = 1, t; = 10, t; = 100, t; = 1000, t; = 10 000, t, = œ, i.e. for: t=0.001s, t = 0.010 s, 
t— 0.100 s, t = 1.000 s, t = 10.000 s, t = œ s, after impulse occurrence the maximum pressure distributions for s, > 0 have the 
following dimensional values, respectively : 


p 0.629228 (66) 


Qo Qo Q 6) (69) Qo 
0.4032. , 0.40325. , 0.4062 , 04342. , 0,588 a , 0,620 2. (67) 
V V V V V V 
CONCLUSIONS 
€ The pressure distribution changes at the instant of impulse  r - dimensional radial coordinate, [m] 
occurrence are caused mainly by the bearing gap height — ri - dimensionless radial coordinate 
changes and viscoelastic oil properties. R - radius of the journal, [m] 
€ The gap height changes during impulsive motion and the Re - Reynolds Number — 
S| - dimensionless coefficient of gap height changes 


viscoelastic oil properties may either increase or decrease 


the pressure distribution and load-carrying capacity of sphe- S ey UL S 


- stress tensor, [Pa] 


rical bearings in contrast to those of the same bearing but gt - Strouhal Number 

free from impulse effects. t - dimensional time, [s] 
€ The pressure distribution changes and load-carrying capa- — t, - dimensionless time 

city values at the instant of impulse occurrence attain about — t, - characteristic value of the dimensional time, [s] 

40 percent of those appearing in the spherical bearing free U - peripheral velocity of spherical journal, [m/s] 

from impulse effects. Vi - dimensional oil velocity component in radial 
€ The pressure distribution changes caused by the viscoela- direction, [m/s] l : l , 

stic oil properties can attain only about 10 percent of those Y= - total dimensionless oil velocity component in radial 

appearing in the instant of impulse occurrence. eoi ; . ; ; 
€ Just after the impulse occurrence the influences of the vi- Vrox - dimensionless oil velocity component in radial 

MA : é direction for Newtonian (classical) unsteady oil flow 
scoelastic oil properties on the pressure and capacity chan- Vix» Vpy...- dimensionless corrections of oil velocity components 
ges quickly tend to zero. This is the moment in which the in radial direction caused by the viscoelastic oil 
largest values of wear may be expected. properties in unsteady flow 

Acknowledgement Vo - dimensional oil velocity component in meridianal 
direction, [m/s] 

This author would like to express his thanks to Polish State ^ voi - total dimensionless oil velocity component in 
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realizing the KBN 8-T 10 B-061-21Grant. Voor - dimensionless oil velocity component in meridianal 

direction for Newtonian (classical) unsteady oil flow 

NOMENCLATURE Voix» Vg2»»-..- dimensionless corrections of oil velocity components 
a - acceleration vector [m/s] in meridianal direction caused by the viscoelastic oil 
A, - strain tensor, [1/s] properties in unsteady flow 
Ay - dimensional tensor, [s] Vo - dimensional oil velocity component in 
Dai - dimensionless origin and end coordinate of circumferential direction, [m/s] 

lubrication surface in meridianal direction Voi - total dimensionless oil velocity component in 
Cy - dimensionless end coordinate of lubrication surface circumferential direction 

in circumferential direction Vox - dimensionless oil velocity component in 
Cii, Cio, Ciz, Cig - integral constant circumferential direction for Newtonian (classical) 

- dimensional eccentricity, [m] unsteady oil flow 

De - Deborah Number Vol» Vo2»»- -.- dimensionless corrections of oil velocity components 
Des = DeStr - Deborah and Strouhal numbers in circumferential direction caused by the viscoelastic 
erf - special integral function oil properties in unsteady flow 
h - gap height, [m] Y,Y,Y,  - dimensionless functions 
hy - dimensional gap height for spherical journal ZZ - dimensionless functions 

and spherical sleeve, [m] a - pseudo viscosity coefficient, [Pas?] 
h; - dimensionless gap height B - pseudo viscosity coefficient, [Pas?] 
Desa - average gap height minimum, [m] € - radial clearance, [m] 
M=Nh,_ - time depended dimensionless function Ae,, A£5, A£3- components of the sleeve centre, [m] 
N - time depended dimensionless function No - dynamic viscosity of the oil, [Pas] 
O - centre of the journal 0 - meridianal direction 
O, - centre of the sleeve (0) - circumferential direction 
p - dimensional pressure, [Pa] x - time depended dimensionless variable 
Di - dimensionless pressure y - dimensionless radial clearance 
Pio - dimensionless pressure for Newtonian (classical) (A) - angular velocity of the journal, [1/s] 

unsteady oil flow 0) - angular velocity of the impulsive changes caused by 


o 


P11» Pi... - dimensionless pressure corrections caused by the perturbations in unsteady conditions, [1/s] 
viscoelastic oil properties in unsteady flow - lubrication surface, [m?] 


re) 
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From 25 to 28 May 2004 the yearly 
21st Symposium on Hydroacoustics 


was held in Jurata, a touristic resort at Hel Peninsula. It was 
organized under the auspices of the Acoustics Committee, 
Polish Academy of Sciences, and the Polish Acoustical So- 
ciety. The symposium was hosted by Naval University of 
Gdynia and Gdansk University of Technology. 


The symposia in question have been aimed at providing 
an opportunity for direct exchange of experience and infor- 
mation among teams dealing with hydroacoustics and rela- 
ted subjects. 


The scope of the 21st Symposium 
covered the following items : 


acoustic wave propagation in sea water 
hydroacoustic noise 

non-linear acoustics in water environment 
ultrasonic transducers 

signal processing 

hydroacoustic devices and systems 

other related problems. 
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Presentation of 31 papers was performed during four ple- 
nary sessions and four topical sessions.The following pa- 
pers were presented during the plenary sessions : 


* Golay's codes sequences in ultrasonography — by A. No- 
wicki, I. Trots, W. Secomski, J. Litniewski (Institute of 
Fundamental Technological Research, Polish Academy 

- of Scieces, Warszawa) 
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Ferromagnetic Oils. Fourth International Congress on Thermal 
Stresses. Proceedings of Thermal Stresses. Osaka, 2001 

11. Wierzcholski K.: Ferromagnetic Turbulent Lubrication for 
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Hydroacoustics 


* Acoustic reconnaisance of fish and environmental back- 
ground in demersal zone in Southern Baltic — by A. Or- 
towski (Sea Fisheries Institute, Gdynia) 

* Underwater ship passport — by I. Gloza (Naval Univer- 
sity of Gdynia) 

* Stability of mechanical and dielectric parameters in pzt 
based ceramics — by J. Ilczuk, J. Bluszcz, R. Zachariasz 
(University of Silesia, Sosnowiec) 

* Directional sonobuoy system for detection of submarines 
by R. Salamon (Gdansk University of Technology). 


Most of the presented papers were prepared by 47 au- 
thors representing 10 Polish universities and scientific cen- 
tres, including the following : of Gdansk University of Tech- 
nology - 8 papers, of institutes of Polish Academy of Scien- 
ces - 6 papers, of Naval University of Gdynia - 5 papers. 
A group of foreign authors consisted of : 3 authors of Se- 
vchenko Research Institute of Applied Physical Problems, 
Minsk, Belarus, 1 author of University of Victoria, Canada, 
1 author of Institute of Applied Physics of Nizhny Novogo- 
rod, and 2 authors of Nizhny Novogorod State University, 
Russia, | author of Institute of Marine Sciences of Mersin, 
Turkey. 
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OPERATION & ECONOMY 


Some aspects of vibration control 


Part I: Active and passive correction 


Krzysztof Kosowski 
Gdansk University of Technology 


ABSTRACT 


The paper presents a general approach to mechanical system modification aimed at con- 
trolling the steady harmonic vibrations by means of passive and active methods. The rela- 
tive decrease of harmonic vibration amplitudes of selected elements of the mechanical 
system has been chosen as a measure of the quality of the introduced modification. The 
proposed theoretical method enables to determine the parameters of the system 5 dynamic 
flexibility matrix, which show the most remarkable effect on the dynamic behaviour of the 
whole system. When active control is considered the method is useful in designing the 


structure and choosing the parameters of the control system. In certain cases of self-excited vibration the 
approach helps examining the elements of the system, most responsible for this kind of excitation. 


Key words : harmonic vibrations, passive control, active control 


PROBLEM DESCRIPTION 


Theoretical investigations into the problem of active con- 
trol of mechanical vibrations have been carried out for many 
years, but real-life mechanical systems making use of this idea 
are still rare. However recently the progress in control methods 
and technology has had an impact on the development of no- 
vel practical solutions, used to reduce vibration of turbomachi- 
nery rotor systems. Application of magnetic and pressurized 
bearings reveals new possibilities for the control of rotor beha- 
viour. With the advent of nanotechnology the practical applica- 
tion of active control methods may increase in the future. The 
application of distributed sensor and actuator systems is tightly 
linked with the theory of multidimensional system control. 

Let us consider a mechanical system with n degrees of free- 
dom, Fig.la. The inertial elements of the system can perform 
translational, bending and torsional movements. The forced har- 
monic vibration of the system is given by the matrix equation : 


J-Q+B-Q+K-Q=F (1) 
where : 

the matrix of the moments of inertia of the system 
the matrix of the damping coefficients of the system 
the matrix of the stiffness coefficients of the system 
the vector ofthe harmonic forces or moments acting upon 
the inertial elements of the system 
the vector of displacement of the inertial elements of the 
system. 


O m"mgEwWwc- 


The matrices J, B, K are of nxn dimension, 
and the vectors F, Q - of nx1 dimension. 
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The forces F acting on the system are assumed to have an 
identical frequency œ, and amplitudes varying within the sys- 
tem to be described by the vector f. They may therefore be 
expressed as a function of time t : 


F=f- e (2) 
where : 
j^ 4-1 - the imaginary unit. 


Fig 1. Schematic diagram of a mechanical system 
and its block diagram describing mechanical vibrations 


The solution Q of equation (1) describing 
the displacement of the elements of the system 
may be expected to have the form : 


Q=q e 


where : 


(3) 


q - the vector of displacement amplitudes 
of the inertial elements of the system. 


After accounting for relations (2) and (3), 
equation (1) takes the form : 


(-w’ J+ joB 4 K)q =f (4) 


All matrices and vectors bearing the subscript i are assu- 
med to refer to the vibrations of the mechanical system for 
a given frequency 0. For the sake of clarity it is helpful to 
introduce the following notation: D; = (K — oJ + jo,B) . The 
matrix Dj; is of nxn dimension. Equation (4) may be rewritten 
to include Dj : 


D; :q; =f; 6 
or equivalently (assuming the equation has a solution) : 
q; - G, ff; (6) 
where : 


the matrix G, = D; stands for the dynamic flexibility 
matrix of the system for the frequency 0o. 


The external forces fz; are assumed to act on selected iner- 
tial elements in a way which can be described by the binary 
matrix Bz;. The matrix qz; in turn represents the vibration am- 
plitudes of the elements whose behaviour has to be controlled. 
The selection of these amplitudes is performed by means of 
the binary matrix Cz;. This is shown (for any vibration frequ- 
ency) in Fig.1b and Fig.1c. 


Active control 


The behaviour of the system shown in Fig.1c is controlled 
by adding a feedback loop. The schematic block diagram of 
the active control of mechanical vibrations is presented in Fig.2. 


Fig 2. Block diagram of active control of mechanical vibrations 


It is here assumed that only the amplitudes qw; of a certain 
number s of inertial elements (selected by the sxn binary ma- 
trix Cw) can be measured. They are treated as an input value 
for the controller. The output vector fw; ofthe controller consi- 
sts of r elements and is a function of the vector qw; of measu- 
red displacement, as well as of the derivative and calculus (of 
any order) of the displacement. The transfer matrix R; of the 
controller action is a rxs matrix with complex elements, which 
fulfils the relation: 


In the case when a PID controller is used in the system, 
the matrix R; takes the following form : 


R; = K, +j (K50; -Kj9 ') 


where : 
Kp, Ky and Kp - rxs matrices of the proportional (P), inte- 
grating (T) and differentiating (D) controller, respectively. 


(8) 


Let us define the rx3s matrix R 
of the controller parameters as : 


R =[K; |K; |Kp] 
The relation between matrices R; and R 
may be written in the form : 


When taking into consideration controller models with dif- 
ferentiation or integration of a higher order it is sufficient to 
extend the parameter matrix R horizontally, and appropriately 
extend the frequency multiplier matrix Q vertically. 

The steering force signal fw; is passed onto the active ele- 
ments of the control feedback loop. The locations of the ele- 
ments upon which they act are given by the nxr matrix Bw. 

The aim of the controller R is to minimize the vibration 
amplitudes of selected elements described by the matrix qzi. 


Passive control 


Passive control is here understood as the modification of 
parameters of the mechanical system. This may be achieved 
by introducing the changes P — [AJ | AB | AK] to certain sys- 
tem parameters (i.e. inertia, damping or stiffness coefficients), 
selected by the binary matrices Bw and Cw. For a given frequ- 
ency 0 the changes introduced to the matrix Dj lead to the 
following changes in equation (5) : 


where : 


P| = PQ; 


(9) 


From equation (9) one obtains : 
qi = D; f, E D; Bw "P, Cw: q; 


By taking into account the notation: G, =D," 
equation (10) may be written in the form : 


q; -G-f; -G-Bw-P. -Cw-q, 


Relation (11) is presented in the form 
of the block diagram shown in Fig 3. 


(10) 


Fig 3. Block diagram of harmonic vibrations 
with correction of system parameters 


The aim of selecting the matrix P; (which describes chan- 
ges of system parameters) is the minimization of the vibration 
amplitudes of chosen elements described by the matrix qzj. It 
should be emphasized that harmonic vibrations with the cor- 
rection of system parameters and the process of active control 
of mechanical vibrations may by represented by the same ge- 
neral form of block diagram (compare Fig.2 and 3). 


Some cases of self-excited vibrations 


In many cases the so-called self-excited vibrations depend 
on the behaviour of the mechanical system itself. Let us consi- 
der, as an example, self excited vibrations of a turbomachinery 
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rotor system due to aerodynamic forces. Rotor-stator eccentri- 
city or rotor-stator misalignment changes the clearance distri- 
bution above the blade shroud (and in glands) which results in 
aerodynamic forces and moments acting on the turbine rotor. 
As a result, self-excited vibrations of the whole rotor system 
may be observed. Due to inaccuracy of manufacture and as- 
sembly, the value and the distribution of the clearance in parti- 
cular seals may differ significantly, and seals of different types 
can be used in the same machine, thus the aerodynamic forces 
acting on a rotor in each stage can vary remarkably. Various 
theoretical models have been elaborated to describe the fluid 
motion in the seals and to determine the aerodynamic forces 
generated in a shroud clearance. Usually the aerodynamic for- 
ces and moments are expressed in the form of rotodynamic 
coefficients [1+ 6, 8, 9, 12, 13] which can be written in the 
form of the following vector equation: 


F=J,Q+B,Q+K,Q 
where : 
F -vector of the components of the aerodyna- 
mic forces and moments 
Q -rotor displacement and rotation in horizon- 
tal and vertical directions 
- matrices with the so called “inertia”, *dam- 
ping" and "stiffness" coefficients of the 
shroud (or gland), respectively. 


(12) 


Ju, By and Ku 


The symbol S is used to denote the matrix 
of rotordynamic coefficients : S=[J, | B, |K,]. 


Because the self-excited vibrations occur at a certain fre- 
quency 0, equation (12) takes a form similar to that of equa- 
tion (4) : 


Bw(Co;J, + jo,B,+K,)Cw-q, =f, (13) 


The matrix Cw selects the displacements responsible for 
aerodynamic excitations, while the matrix Bw determines the 
places where these forces are applied. 

The equation (13) can also be written in the form : 


where : 
S, 2-0;J, + jo;B, +K, =S-Q, 


When taking into consideration only the aerodynamic exci- 
tation forces, the behaviour of the rotor system can be repre- 
sented by means of the block diagram shown in Fig.4. 


Fig 4. Block diagram of harmonic self-excited vibrations 


The comparison of Fig.2 , 3 and 4 leads to the conclusion 
that harmonic vibrations with the correction of system para- 
meters, the process of active control of mechanical vibrations 
as well as some cases of self-excited vibrations can by repre- 
sented by one and the same general form of block diagram. 
Thus all the three cases may be generalized to the form presen- 
ted in Fig.5. The matrix U; can represent any of the matrices 
R;, Pj, Sj, depending on the context. All further considerations 
are conducted with the use of this general form, and the atten- 
tion is focused on determining the value of the coefficients of 
the matrix U for which the amplitudes qz of chosen elements 
of the system (described by the dynamic flexibility matrix G) 
will be minimum. 
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Fig 5. Block diagram of the system performing mechanical vibrations 


When taking into consideration vibrations of a mechanical 
system with k different values of frequency œ it is possible to 
generalize the model given in Fig.5 to the form shown in Fig.6. 
This is performed by introducing generalized matrices G, Q, 
U, U, Bw, Cw, Bz, Cz, and vectors f, fz, fw, q, qz, qw (Fig.6), 
which correspond, respectively, to the matrices Gi, O;, U, Ui, 
Bw, Cw, Bz;, Cz;, and vectors f, fz;, fwi, qi, qzi, qw; in the case 
of vibrations with a single frequency œ; (Fig. 5). 


Fig 6. Block diagram of active control of mechanical vibrations 


All of the vectors f, fz, fw, q, qz, qw in the generalized 
scheme are column vectors composed of the corresponding 
vectors for a single frequency. In other words, these vectors 
may be written in the form : 


f fz fw 
f=| ?| fz— .°| fw=| .° 
f, fz, fw, 
qı qZ; qW; 
1q? || dZ5 ELE 
q= E qz — : qw — 5 
qk qZ, qw | 


The matrices G, Q, U, U, Bw, Cw, Bz, Cz 
are described by the following formulas : 


G 0 -.. 0 Qi 0 0 
"E Gyr 0 |o- 0 Q, - 0 
0 0 G, 0 0 ux 
U 0 0 U 0 — 0 
e Un ofp jeune 
0 0 U 0 0 : Uk 
Bw 0 0 Cw 0 = 0 
Bw- 0 Bw 0 Cw-| 9 Cw - 0 
0 0 Bw 0 0 -> Cw 
Bz, 0 0 Cz, 0 0 
ul) mo jet emen 
0 0 Bz, 0 0 Cz, 


It is useful to note that all the generalized matrices may be 
described as a linear combination of the corresponding matri- 
ces for particular frequencies and certain other fixed matrices. 
In particular, for the matrix U it can be written : 


k 
Ü- Y e xr)-U-e?xr,)] 


i=l 


(14) 


where : 
e; - the i-th versor of dimension k, the symbol : x stands for 
the carthesian product of matrices. 


METHOD OF RESPONSE CIRCLES 


The following set of equations can be written 
for the system presented in Fig.6 : 


q=G-f (15) 
f = Bz -fz + Bw fw (16) 
qz = Čz -q (17) 
qw =Cw-q (18) 
fw =U-qw (19) 
U-Ü.o (20) 
By combining equations (14) + (22) 
the following relation can be derived : 
qz= Cz(I,- G-iwY: ond | Q. Cw)" G-fz 
izi |: U - (e; x Ly) (21) 


In the case without any feedback (U=0) the vector qzo 
of amplitudes of selected elements may be written in the form : 


qzo = Cz: G: fz (22) 


By using equations (15) + (21) it is possible to investigate 
how the real or imaginary part u ofa particular coefficient of the 
matrix U influences the amplitudes in vector qz. The vibra- 
tions of particular elements represented by the matrix qz can 
play a different role in the dynamic behaviour of the mechani- 
cal system. Therefore in some cases it is useful to use a weigh- 
ted sum of amplitudes as a measure of the vibration level : 


o, = 2 ali] -qzli] 


where : 
a[1] - the weight coefficient 
corresponding to amplitude qz[i]. 


23) 


For the case of no feedback it can be similarly written : 
6,9 = Y ali]-qz,[i] 
i 


After some transformation it is possible to show that the 
ratio G of the indexes o; and 079 may be written in the general 
form : 


(24) 


,(u) K u(a; + jb;) 
Qu) = S29 214 y —RU n - 
O, (U) ii 1— u(c; + jd;) 
where : 
aj, bi, cj, dj, (1 <i € k) are real numbers. 


The module |C| shows the effect of the feedback loop on the 
weighted vibration amplitude of the elements of qz. In gene- 
ral, when V... [d; 0^ (a; € 0v b; #0)] equation (25) repre- 
sents a closed smooth curve which is described by the end of 
the vector G in the complex coordinate system (r , i) when u 


Q5) 


varies from - œ to + co, An example of such curve is presented 
for k = 4 in Fig.7. 


Fig 7. Graphical interpretation of equation (25) 
with exemplary coefficients for k different frequencies (k = 4 ) 


When taking into consideration only a particular 
frequency w, equation (25) may be written in the form : 
u(a + jb) 


Gu) =1+ 
1 - u(c 4 jd) 
where : 
8, b, c, d, are real numbers. 


(26) 


In the usual case when d z 0 A (a z 0 v b z 0), equation (26) 
represents a circle which is described by the end of the vector G 
in the complex coordinate system (r , i) when u varies from - eo 
to + œ [7] : 


Fx] ia 


A graphical interpretation 
of equations (26 , 27) is shown in Fig.8. 


a) b) 


a^ b? 


eC 


i 


Fig 8. Graphical interpretation of equation (27) 


The diameter of the closed curve (in particular — a circle) 
described by the vector G enables to assess the influence of the 
parameter u on the weighted amplitude of vibrations. The large 
diameter of the curve (the circle described by the vector C) in 
Fig.8a) may be interpreted as a significant influence of the pa- 
rameter u, while the small diameter (the circle determined by 
Co in Fig.8a) — as an insignificant influence of the parameter u. 
However in practice the parameter u has a reasonably limited 
range of values. Thus only a part of the curve drawn by the 
vector G can be applied in practice (for example only the mar- 
ked part of the circle shown in Fig. 8b). Moreover, only the 
values of parameter u, for which the module |C] is lesser than 1, 
result in the decrease of the weighted amplitude level. In this 
way it is possible to estimate the effect of all parameters u on 
the vibrations and thus to choose : 


* the controller structure and its parameters for the best acti- 
ve control of mechanical vibrations 
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* the changes of system parameters leading to the most ef- 
fective reduction of vibrations 

* the shrouds and the glands which play the most important 
part in generating forces and moments in the case when 
rotor self-excited vibrations of aerodynamic type are con- 
sidered. 


EXAMPLES 
A ship propulsion system 


Ship propulsion systems equipped with flexible couplings 
are very sensitive to disturbances caused by unsteady engine 
operation. The disturbances have the form of shaft torque pe- 
riodical changes which lead to torsional vibrations of the who- 
le propulsion system. In some cases resonance vibrations re- 
sulting in damages to flexible couplings, were observed. This 
situation very often occurs when the engine works with one 
misfiring cylinder. In Fig.9a an example ship propulsion sys- 
tem is presented. 


a) 


GENERATOR 


MAIN ENGINE 


c) 
[ H Ho 
Ho 

d) 
[ H HC 
Ho 


Fig 9. Variants of additional control systems and their response circles 


It consists of a main medium-speed diesel engine which — 
— through a main coupling and a mechanical gear — drives 
a ship propeller, an electric generator and two hydraulic pumps. 
All the couplings are flexible. The linear model of this system 
and the analysis of its behaviour was elaborated [11] by using 
the engine producer’s data and results of some additional inve- 
stigations. The reduction of the main coupling torsional vibra- 
tion was performed by modifying the main engine governor. 
Three system variants were considered for the following diffe- 
rent correction input signals : 


> angular velocity of the generator (Fig.9b) 
> angular velocity of the main coupling (before the gear) 
(Fig.9c) 
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» angular velocity of the propeller shaft (measured directly 
after the gear) (Fig.9d). 


In all three cases the presented method of response circles 
was applied to the analysis of the structure and parameters of 
the additional correction controller. The reduction ratio G of the 
main coupling vibration amplitude for the case of a proportio- 
nal controller, 270 rpm shaft speed and 14 Hz fundamental har- 
monic frequency of forced vibrations, is shown in Fig.9b + 9d 
for the above mentioned system variants, respectively. 

From Fig.9 it is evident that the controller using the gene- 
rator's angular velocity as its correction signal offers the lar- 
gest possibilities of reducing torsional vibration amplitude in 
the main coupling. 


Turbine rotor self-excited vibrations 


The method of response circles was used to select the seals 
which play the most important part in generating aerodynamic 
forces leading to self-excited vibrations of the rotor system. 
The forces were described by means of the rotordynamic coef- 
ficient matrix S=[J, | B, | K, ], and calculated from relation 
(13). The relevant schematic diagram is shown in Fig.10. 


E 


au} | m 
P 
4 


S 


Fig 10. Schematic diagram for rotodynamic coefficient analysis 


In the case of a double-cylinder medium-power steam tur- 
bine the performed analysis proved that the seals of the first 
stages of the HP cylinder had the greatest influence on self- 
-excited vibrations of the aerodynamic type. It was enough to 
change the seals of the shrouds and shaft in first four turbine 
stages to achieve the desired effect of vibration reduction. 

A currently conducted work is concentrated on active con- 
trol of rotor vibrations of a steam turbine by means of pressu- 
rized bearings. The method of response circles is used to de- 
tect the bearings which have the greatest influence on active 
control. Results of the work in question will be presented in 
a separate paper in due course. 


CONCLUSIONS 


O An analytical method for the investigation of linear mecha- 
nical systems performing harmonic motion was presented. 
This approach was successfully applied for the analysis and 
improvement of the dynamic behaviour of a ship propul- 
sion system. 

O The proposed theoretical method makes it possible to de- 
termine the parameters of the system’s dynamic flexibility 
matrix which show the most remarkable effect on the dy- 
namic behaviour of the whole system. 

O When active control is considered the method is useful in 
the designing of the structure and choice of parameters of 
the control system. 

O Incertain cases of self-excited vibrations the approach helps 
examining the elements of the system which are most res- 
ponsible for this kind of excitation. 


In Part II of the paper (to be published) a theoretical me- 
thod of optimum vibration control by active means is descri- 
bed and illustrated by some examples, including an approach 
to pressurized bearing application. 
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NOMENCLATURE 

B - matrix of damping coefficients 
C - selection matrix 

Bw, Bz, Cw, Cz - binary matrices 

D - transfer matrix 

f - vector of force amplitudes 

fw - controller output vector 


- vector of external forces 

- vector of harmonic forces (or moments) 

- dynamic flexibility matrix 

- unitary matrix 

- imaginary unit 

- matrix of inertia moments 

dimensions of matrices and vectors 

- matrix of stiffness coefficients 

- matrix of changes of system parameters 

- vector of displacement amplitudes 

vector of measured amplitudes 

vector of amplitudes of controlled elements 
- vector of displacements 

- controller matrix 

- matrix of rotordynamic coefficients J, , B, , K, 
- time 

- general symbol for the matrices P, R, S 


BP 

I 

un 
1 


N 
1 


- weighted sum of amplitudes 
- frequency 
- frequency multiplier matrix 


psa d^^xuOoqSomvRreccomug 


Indices: P - proportional controller 
I -integrating controller 
D - differentiating controller 
u -ofrotordynamic coefficients of turbine seals 
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On possible supplying ship diesel engines with 
alternative fuels (mixtures of diesel oils and 
vegetable oils or their esters) 


Preliminary report 


Jacek Krzyzanowski 
Kazimierz Witkowski 
Gdynia Maritime University 


ABSTRACT 


The paper presents introduction to the research on possible supplying ship diesel engines with mixtures of 

diesel oils and vegetable oils or their esters with accounting for ecological aspects, i.e. exahust gas purity. 

Characterisitics of vegetable oils and their esters are compared with those of diesel oils; some consequen- 

ces of their application to diesel engines, mainly for their working process and exhaust gas content, are 

indicated. Also, influence of combusting their mixtures with diesel oils are discussed in the same context. 
Scope of the planned research project is shortly presented. 


Key words : ship diesel engines, alternative fuel oils, ecology 


INTRODUCTION 


Contemporary main diesel engines of sea-going ships are 
commonly supplied with heavy oil fuels. This very often con- 
cerns also auxiliary engines, especially electric generating sets. 

However on many ships the electric generating sets are still 
fed with marine diesel oils (MDO). Also, most of diesel engi- 
nes installed on small ships are run on MDO. 

Permanently increasing demand of diesel oils (DO), incre- 
ase of their prices, increasing ecological requirements make that 
more and more attention is paid to the alternative fuels called 
also substitute, renewable or unconventional. Another reason 
of the growing interest to the fuels is the increasing probability 
of dropping worldwide output of crude oil due to different cau- 
ses. Moreover the problem of the hazard to natural environment, 
associated with mining, transport, processing and combusting 
the oil products is today brought up more and more strongly. 

All energy sources other than crude oil products may be 
deemed unconventional ones [1]. As it results from the dia- 
gram presented in Fig.1. there is a wide range of possible me- 
dia for supplying diesel engines. However today most of the 
unconventional fuels are not used in practice of operation of 
diesel engines, including ship engines. 

In the research project planned by these authors the supply- 
ing of a ship diesel engine only with a mixture of marine diesel 
oil (MDO) and rape-feed oil methyl ester (RME) has been 
accounted for. 


Application of vegetable oils 
for supplying diesel engines 


In the northern zone of moderate climate rape oil, flaxseed 
oil and corn oil are mostly taken into account. In other climatic 
zones it can be soybean, sunflower, palm, cotton, sesame, pea- 
nut or coconut oil. 
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Unconventional fuels 
for supplying diesel engines 


Compressed Liquefied 
Natural Gas 
(earth) gas Liquefied 
Digester | TING) 
(bio)gas 


and mineral oils 
Town i 


(iluminating) gas l Fatty Acid 
; Methyl 
Hydrogen H Esters (FAME) 


Propane | |  ., " - 
mixture Diethyl ether 


Fig. 1. Scope of unconventional fuels for supplying diesel engines [1] 


Vegetable oils are esters of glycerol and fatty acids contai- 
ning from 14 to 22 carbon atoms [5]. In Tab.1 [1] some proper- 
ties of diesel oils and selected vagetable oils are compared. 

It can be observed that some parameters of vegetable oils, 
their esters and diesel oils are of similar values, respectively. 
However their density, kinematic viscosity and flow tempera- 
ture values much differ from those of diesel oils. 

As far as rape oils are concerned (most interesting in case 
of Poland) their density and viscosity is distinctly higher, which 
can make supplying diesel engines with them difficult, how- 
ever their positive features are: practically no sulphur content 
and bio-degradation ability. Their calorific value is smaller than 


Tab. 1. Comparison of some properties of diesel oils, vegetable oils and methyl esters of higher fatty acids of rape oil (RME) 
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that of diesel oils, due to different chemical content. Results of 
research on application of vegetable oils for supplying diesel 
engines [2, 3, 4, 9, 10] show worse cylinder filling, worse spray- 
ing, greater lengths of injected oil jets, associated with their 
large viscosity and density. Engines supplied with rape oil ope- 
rate with a lower total efficiency mainly due to its lower calo- 
rific value in comparison with that of diesel oils, as well as due 
to worse spraying process resulting in a longer combustion. 
There is no unambiguous results of research on exhaust gas 
toxicity. Most obtained data indicate an improvement of exhaust 
gas purity relative to its content in the case of supplying engi- 
nes with diesel oil. However the phenomena associated with 
supplying diesel engines with rape oil are disturbing, namely : 


* often occurence of clogging sprayer nozzles 
* great susceptibility to forming carbon deposits on piston 
heads, ring grooves, valves and valve seats 


According to [1,10] 20% addition of diesel oil to rape oil 
made its viscosity dropping by 30% as well as ignition lag pe- 
riod lowering; engine starting features appeared improved. In 
common operational applications small additions of rape oil or 
RME to diesel oil (e.g. 5% to 20%) are usually reported. Tests 
on mechanical vehicles running on a 20% RO / 80% DO mix- 
ture did not reveal any detrimental consequences [10]. 

Therefore these authors have decided to carry out investi- 
gations on diesel engines supplied with MDO/RME mixtures 
initially containing no more than 10% of RME. 

The basic properties of the MDO and RME and their mixtures 
selected for the experiments are given in Tab.2. 

Preparation of such mixture is easy as it does not reveal 
a tendency to separation. However such mixtures could be less 
ressistant to ageing process hence their should be consumed in 
a short time. 


Tab. 2. The basic properties of the MDO and RME and their mixtures selected for the experimental tests 


Viscosit 


50°C 
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Rape oil methyl 
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95% MDO + 
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* troubles with starting engines at low ambient temperatures 
* seizing precise pairs of injection pumps. 


For these reasons it seems not advisable to commonly use 
vegetable oils as fuels for diesel engines. However much bet- 
ter effects can be obtained by applying chemically modified 
oils, i.e. methyl esters of higher fatty acids (FAME), see Fig.1. 
In Poland and many other European countries rape oil and 
methyl alcohol are usually applied to produce the esters called 
rape-feed oil methyl esters (RME). 


Mixtures of diesel oils 
and vegetable oils or their esters 


Due to substantial difficulties in applying only vegetable 
oils as well as due to limitation in using esters for running die- 
sel engines, an alternative is to use mixtures of diesel oils and 
vegetable oils or diesel oils and vegetable oil esters. Properties 
of such mixtures depend on properties of their components and 
their content in a mixture. In this way it is expected to decrease 
density and viscosity of a mixture relative to those of a given 
vegetable oil or ester. 


SCOPE OF THE PLANNED 
RESEARCH PROJECT 


The planned research project is aimed at revealing : 


+ the influence of using MDO/RME mixture to running 
a ship diesel engine on its operation 

+ other possible problems arising from using such fuel mix- 
ture, e.g. its durability, storage conditions, transport, filtra- 
tion processes etc. 


The investigations will be carried out with the use of the 
L22 diesel engine installed in the laboratory of ship combu- 
stion engines, Gdynia Maritime University. The laboratory 
engine and test stand to be used were described in [6, 7, 8]. 

The planned tests will have a comparative character.The 
first series of all the tests will be devoted to the investigations 
of the engine fed with pure marine diesel oil, results of which 
will be taken as model ones. The consecutive test series will be 
carried out with MDO/RME mixtures of various content, and 
the controlled parameters such as engine speed or load will be 
changed in the same way as in the model series. Both absolute 
and relative values of the parameters will be analyzed. 
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The largest group of ship diesel engines working on diesel 
oil are those driving electric generating sets, which, due to the- 
ir function, are of constant speed and highly changeable load 
(torque or fuel oil charge). Therefore the first tests will be per- 
formed for such loading mode of the engine. At different but 
constant engine speed values the engine load will be changed 
within the range from 20% to 80% of the rated torque, during 
which all important parameters of the engine working process 
as well as indication diagrams will be recorded. During the 
tests special attention will be paid to exhaust gas content both 
from the point of view of engine working process course and 
noxious component content. 


In the first phase of the research the use of two fuels (of 5% 
and 10% of RME content in MDO/RME mixture) is assumed 
because in the case no changes in the engine’s construction or 
regulation system are necessary. 


The next phase of the research should be focused on durabi- 
lity (reliability) of the engine installed on a ship. However deci- 
sion on initiating such tests will be taken on the basis of analy- 
sis of the results obtained in the first phase of the research. 
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~~~ C onference T 
EXPLO-SHIP 2004 


On 31 May and 1 June 2004 3rd International 
Scientific Technical Conference on : 


Operational Problems of Floating Objects 
and Port Facilities 


was held in Świnoujście, a Baltic coast 
port and health resort. 


It was organized by Sea Navigation Institute, 
Maritime University of Szczecin. 


During its scientific part 124 papers including 
the following 5 plenary ones, were presented : 


* Coordinate systems applied to navigation over limited 
water regions — by S. Gucma (Maritime University of 
Szczecin 

* Some proposals on an unified maritime safety and se- 
curity system — by Z. Kopacz, W. Morgas, J. Urbanski 
(Naval University of Gdynia) 

* Operational reliability of transport systems — by J. Jaz- 
wiiski, Z. Smalko, J. Żurek (Air Force Institute of 
Technology, Warszawa) 

* Prospects of application of natural gas to ship engines 
by St. Zmudzki (Technical University of Szczecin) 

* Contemporary systems for determining search areas 
at sea — by Z. Burciu, J. Soliwoda, S. Ukleja (Gdynia 
Maritime University) 


The remaining 119 papers were presented 
during 8 topical sessions on : 


Modelling and simulation 

in navigation (15 papers) 

Operation of ship (13 papers) 
Navigational systems (16 papers) 
Safety at sea (12 papers) 

Ship combustion engines (14 papers) 
Ship systems (17 papers) 

State identification (16 papers) 
Operation of ship systems (16 papers). 


+++++4+4+ + 


Several authors of the papers came from Russia and 
Ukraine. An attractive end of the Conference was the tou- 
ristic voyage of its participants to Copenhagen, Denmark, 
onboard a ferryship. 


C nfi —N 


PIAP - Automation 2004 


On 24-26 March 2004 
the Technical Scientific Conference on : 


Automation - novelties and prospects 


organized by the Industrial Institute for Automation 
and Measurements, was held in Warszawa. 


It was already the 8th - since 1997 - meeting of experts 
from scientific, research and development centres, as well 
as industrial enterprises, which was a good opportunity to 
present achievements and to exchange experience in the area 
of practical applications of means for automation and robo- 
tics as well as of measurement instruments and systems ap- 
plicable to various fields of engineering. 


The Conference program comprised 4 plenary papers 
and 61 papers to be presented during 
six topical sessions as follows : 


I - Automation, robotics, monitoring (34 papers) 

II - Software, equipment and applications of mobile robots 
(7 papers) 

III - Methods of design and integration of systems 
(4 papers) 

IV - Devices for automation and robotics (9 papers) 

V - Measurement instruments and systems (5 papers) 

VI - Industrial network communication systems (2 papers). 


The authors of the papers represented over 20 universi- 
ties, scientific institutes and research centres, among which 
specialists from Czech Republic and Germany were also pre- 
sent. The greatest number (14) of the papers were prepared 
by the authors from Warsaw University of Technology, and 
9 and 8 - by the authors from the Industrial Institute for Au- 
tomation and Measurement, Warszawa, and Silesian Uni- 
versity of Technology, Gliwice, respectively. 

The following papers, directly or indirectly dealing with 
maritime engineering and economy, were presented in the 
plenary session : 


* Methods and systems for monitoring industrial proces- 
ses — by J. M. Koscielny (Warsaw University of Techno- 


logy) 


* Synergy in mechatronics — by Z. Gosiewski (Military 
University of Technology, Warszawa) 


Session I : 


* PID or fuzzy logic control of water level in steam 
boiler? — by P. Biały and S. Skoczowski (Technical Uni- 
versity of Szczecin) 

* Dynamic decoupling of right-hand-side convertible dy- 
namic systems — by P. Dworak and S. Banka (Technical 
University of Szczecin) 

* Steadfast MFC-PID temperature controller and its re- 

alization by means of PLC programmer — by K. Pie- 

trusewicz and S. Skoczowski (Technical University of 

Szczecin 


* Flow control of pallets applied in a research flexible 
manufacturing system — by J. Honczarenko and M. Sos- 
nowski (Technical University of Szczecin) 


Session II : 


* Diagnostics of failure states of an underwater vehicle in 
service conditions — by T. Leszczyński (Polish Naval Uni- 
versity, Gdynia) 

* Automatic control of immersion depth of an underwa- 

ter vehicle — by P. Szymak (Polish Naval University, 

Gdynia) 


Session III : 


* Aiding of design of ship power plant automation systems 
by using of CBR method — by M. Meler-Kapcia, 
Z. Kowalski and S. Zielinski (Gdansk University of Tech- 
nology) 


Session IV : 


* Computer testing of automatic synchronizers for elec- 
tric generators — by A. Grono, G. Redlarski and J. Za- 
walicz (Gdansk University of Technology) 

* Conical magnetic bearings — by Z. Gosiewski, K.Fal- 

kowski (Military University of Technology, Warszawa) 

and L. Matuszewski (Gdansk University of Technology) 


Session VI : 


* Application of USB buses to microprocessor automation 
systems — by M. Porzezinski and M. Drewka (Gdansk 
Universty of Technology). 


The Conference participants were also given the oppor- 
tunity to visit 10th International Fair for Measurement and 
Control and 5th International Fair for Pumps and Industrial 
Fittings. 


FOREIGN 


AIWARM 2004 


On 26-27 August 2004 
the Asian International Workshop on : 


Advance Reliability Modelling 


had place in Hiroshima City, Japan. During 20 topical ses- 
sions 79 papers were presented. Their authors came from 
scientific centres of Canada, China, India, Italy, Japan, Ko- 
rea, Kuwait, New Zealand, Russia, Singapore, South Afri- 
ca, Sweden, Taiwan, Thailand and USA. 

Prof. K. Kolowrocki, a representative of Gdynia Mariti- 
me University (Poland) was among them, who presented 
the ordered paper on : 


Reliability and risk evaluation of large systems. 
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FOREIGN 


Science - Practice- Education 


Under this heading scientific cooperation among me- 
chanical faculties of Polish and German universities is 
developed. Its starting point was a symposium held in 
Bremen in 1989, in which representatives of local univer- 
sities and scientific workers from Mechanical Faculty, 
Gdansk University of Technology, took part. 

With time, the yearly meetings during which experience 
and scientific information has been exchanged on the ba- 
sis of presented papers, has stirred up greater and greater 
interest. Both number of participants and group of co-or- 
ganizing universities has been growing and program of 
the symposia enriched by visits to industrial enterprises 
and participation in local jubilee ceremonies. 


On 5+7 May 2004 14th Symposium of the kind, 
organized by a German Technical University, 
was held in Stralsund. 


From Polish side 8 persons from Gdansk University 
of Technology, 4 - from Technical University of Szczecin, 
3 - from the State Higher Engineering School of Elblag, 
and | - from the Fluid-Flow Machinery Institute, Polish 
Academy of Sciences, Gdansk, took part in the Sympo- 
sium. 


The following papers were presented 
by Polish participants : 


from Gdansk University of Technology 


9* Investigations of wall influence on single bubble mo- 
vement — by D. Mikielewicz and J.Wajs 

* Hydrodynamic Tunnel - the test stand for investiga- 
tion of characteristics of profiles and cascades 
by J. Iwan, B. Lednicka, K. Zochowski 

* New types of hydraulic pumps and motors prepared 
by the Chair of Hydraulics and Pneumatics, Gdansk 
University of Technology, for practical applications 
by A. Balawender, L. Osiecki 


from Technical University of Szczecin 


# Carbon nanoforms - a new material in material science 
by R.J.Kalenczuk, E. Borowiak-Palen 

# Evaluation of geothermal heat utilization possibilities 
in heat plant cooperating with heat receivers connec- 
ted in parallel — by W. Nowak, K. Zwarycz 


from the State Higher Engineering School of Elblag 


* Dynamic analysis of high-power turboset — by M. Kah- 
sin, H. Olszewski, Z.Walczyk. 


ASME - TURBO EXPO 2004 


From 14 to 17 June 2004 a global gas turbine event for 
education, technology and networking had place in Vien- 
na. It was the international Technical Congress of Ameri- 
can Society of Mechanical Engineering (ASME) under the 
heading : 


Power for Land, Sea and Air 


During 214 sessions were presented papers prepared 
by almost 2000 specialists representing universities, scien- 
tific institutes, scientific research centres and industrial 
enterprises of many countries of the world. In the event 
also Polish scientific workers took part. They presented 9 
papers as follows : 


* Mathematical model of Solid Oxide Fuel Cell (SOFC) 
for power plant simulations — by A. Miller, J. Milewski 
(Warsaw University of Technology) 

* Optimization of rotor critical speeds by change of 
features of machine's bearings — by J. Rybczyński 
(Institute of Fluid-Flow Machinery, Gdansk, Polish 
Academy of Sciences) 

% Experimental analysis and prediction of wake-indu- 
ced transition in turbomachinery — by W. Piotrowski, 
S. Drobniak, W. Elsner (Technical University of Cze- 
stochowa) and S. Vilmin (Numeca International, Bel- 
gium) 


* Natural frequencies and modes shapes of two tuned 
and mistuned bladed discs on the shaft — by R. Rzad- 
kowski, M. Drewczynski (Institute of Fluid-Flow 
Machinery, Gdansk, Polish Academy of Sciences) 

* Unsteady load of partial admission control stage ro- 
tor of a large power steam turbine — by R. Rzadkow- 
ski, P. Lampart, M. Szymaniak (Institute of Fluid-Flow 
Machinery, Gdansk, Polish Academy of Sciences) 

* Tip leakage/main flow interactions in multi-stage HP 
turbines with short-height blading — by M. Szymaniak, 
P. Lampart (Institute of Fluid-Flow Machinery, Gdansk, 
Polish Academy of Sciences) and S. Yershov, A. Ru- 
sanov (Institute for Mechanical Engineering Problems, 
Ukraine) 

* Application of ANN for diagnostics of the geometry 
deteriorations of the power system apparatuses 
by J. Gtuch (Gdansk University of Technology) and 
J. A. Krzyzanowski (Institute of Fluid-Flow Machine- 
ry, Gdansk, Polish Academy of Sciences) 

* An impact of environmental disturbances on combi- 
ned cycle power plant control — by Z. Domachowski, 
M. Dzida (Gdafisk University of Technology) 

* Expansion line modelling and strength diagnostics 

of internally cooled gas turbines — by W. Kosman, 

T. Chmielniak (Silesian University of Technology, Gli- 

wice). 
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